CHEMICAL ENGINEERING SCIENCE 
GENIE CHIMIQUE 


1959 








Grossen, Einheiten und Einheitensysteme 


U. Gricunus 


Farbenfabriken Bayer A. G., Leverkusen, W. Germany 
( Received 11 December 1958) 


Zusammenfassung —Dic Beschliisse der 10 Generalkonferenz fiir Mass und Gewicht (1954) 
sowie die Empfehlungen der © International Organisation for Standardization” (ISO) und anderer 
internationaler Fachorganisationen zeigen das Bestreben, zu grésserer Einheitlichkeit in dem 
Gebrauch von physikialschen Gréssen und ihren Einheiten zu gelangen. In der Bundesrepublik 
Deutschland hat der Verein Deutscher Ingeneiure kirzlich Empfehlungen ausgesprochen, die dem 
uleichen Ziel dienen sollen und tiber die im folgenden berichtet wird. 

Abstract—The resolutions of the 10th General Conference on Weights and Measures (1954) and 
the recommendations of the International Organisation for Standardization (ISO) and of other 
international organisations indicate the desire to achieve greater uniformity in the use of physical 
quantities and their units. In the German Federal Republic the Verein Deutscher Ingenicure has 
recently made recommendations intended to serve the same purpose, which are reported cn in 


the following. 


Mesures (1954), les 
OIN) et d'autres 


organisations internationales manifestent le desir dachever une plus grande uniformitée dans 


Résumé Les résolutions de la ioe Conférence Générale des Poids et 


recommendations de FOrganisation Internationale pour la Normalisation (ISO 


remploi de grandeurs physiques et de leurs unites. Pour la République Fédérale Allemande. 


auteur indique les recommandations rédigees par le Verein Deutscher Ingenicure. 


der Grad Kelvin (K) als Einheit der thermo- 
dynamischen Temperatur, 
die Candela (cd) als Einheit der Lichtstdrke. 


1. EINLEIrTUuNG 
Mir grosser Aufmersamkeit wird in allen Kultur- 
staaten die neuere Entwicklung auf dem Gebiet 


der Einheiten und Einheitensysteme verfolgt, die : :; oe " 
ae —*< = Das zu diesen Grundeinheiten koharente Einhei- 
durch die Beschliisse der 10. Generalkonferenz ti <a 

tensystem wird im folgenden mit “Internationales 


fiir Mass und Gewicht vom Jahre 1954 ausgelést ; = ; ; 
. ‘“ Linheitensystem bezeichnet. Haufig  spricht 


Diese Konferenz setzte 6 Grundeinheiten — « 
: man auch abkiirzend vom MKS-System oder vom 


wurde. 


fest und empfahl die allgemeine Verwendung des 
durch Kin- 
heitensystems in den Landern der Meterkon- 


diese Grundeinheiten bestimmten 
vention. 


Diese 6 Grundeinheiten sind: 


das Meter (m) als Einheit der Lange, 

das Kilogramm (kg) als Einheit der Masse, 

die Sekunde (s) als Einheit der Zeit, 

das Ampere (A) als Einheit der elektrischen 
Stromstarke, 


MKSA-System. 

Andere — internationale Fachorganisationen 
schlossen sich dem Vorgehen der Generalkonferenz 
* International 
Standardization ” (ISO), die “ International 
Union for Pure and Applied Physics (IUPAP), 
die * International Union for Pure and Applied 
Chemistry ” (IUPAC) 
Electrical Commission” (TEC). 


an, so die Organization for 


* International 
In vielen Doku- 
menten dieser Organisationen, die auch Einheiten 


und die 
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anderer Systeme oder systemfreie Einheiten 
enthalten, ist das Internationale Einheitensystem 
deutlich hervorgehoben und empfohlen. 

Angesichts dieser Entwicklung mussten sich 
auch die deutschen Ingenieure die Frage vorlegen, 
ob ein weiteres Festhalten am metrischen tech- 
nischen Einheitensystem verniinftig sei. Bei der 
zunehmenden Verflechtung der Ingenicurwissen- 
schaften mit Physik, Chemie und physikalischer 
Chemie, wie sie gerade fiir die Verfahrenstechnik 
(Chemical Engineering) typisch ist, erwies sich 
das technische Einheitensystem als binderlich fir 
eine gegenseitige Verstandigung. 

Das metrische technische Einheitensystem 
beruht (im Bereich der Mechanik) auf 
Grundeinheiten Meter, Kilogrammkraft und Sek- 
Inter- 


unde. Es unterscheidet sich also vom 
nationalen Einheitensystem oder vom (hier nicht 
weiter behandelten) CGS-System nicht nur durch 
andere Grundeinheiten, sondern auch dadurch, 
dass eine diesser Grundeinheiten zu einer anderen 
Gréssenart gehérte. Besonders verwirrend wirkte 
Kilogrammkraft durchweg mit 


es, dass das 


Kilogramm bezeichnet wurde, obwohl bereits die 


1. Generalkonferenz (1889) das Kilogramm 
eindeutig zur Einheit der Masse erklirte und 
seither diese Festlegung wiederholt bestatigt wurde. 

Die oben genannte Frage, ob das metrische 
technische Einheitensystem im deutschen tech- 
nischen Schrifttum beibehalten werden sollte, 
wurde durch den Wissenschaftlichen Beirat des 
Vereins Deutscher Ingenienure (VDI) diskutiert 
und nach langen Beratungen schliesslich verneint. 
Die Beratungen fanden dadurch ihren Abschluss, 
dass Empfehlungen, die eine hierzu eingesetzte 
Kommission augesarbeitet hatte, zusammen mit 
einem Kommentar von FLEGLER [1] verOffentlicht 
wurden. 

Im folgenden sind diese Empfehlungen im 
Wortlaut wiedergegeben. Einzelne Punkte, soweit 
sie insbesondere fiir die Leser dieser Zeitschrift 
von Interesse sein kénnen, werden ausfihrlicher 
besprochen. 


2. EMPFEHLUNGEN VDI 
1. Die in Z. Ver. dtsch. Ing. 1950 92 161 
veroffentlichten Stellungnahmen Wissen- 
schaftlichen Beirates sind aufzugeben. 


DES 


des 


2. Unter Kilogramm (kg) ist in Zukunft nur 
noch das Massenkilogramm zu verstehen. Der 
Krafteinheit im Technischen Mass-System_ ist 
eine andere Bezeichnung als Kilogramm, z.B. die 
Bezeichnung Kilopond (kp), zu geben. 

3. Das Internationale Einheitensystem mit 
den sechs Grundeinheiten Meter, Kilogramm, 
Sekunde, Ampere, Grad Kelvin, Candela, ist zu 
bevorzugen. Fir dieses Einheitensystem gilt : 


3.1 Die Einheit der Kraft ist das Newton (N), 
IN = Ilkgms™. 


Das Kilopond ist damit definiert durch die 
Gleichung 

Lkp = 9,80665 N, 

Die Einheit (Arbeit, 
Warmemenge usw.) ist das Joule (J), 


1J =1Nm. 


3.2 der Energie 


Als weitere Energie-EKinheiten kénnen ver- 
wendet werden : die Kilowattstunde (kWh), 


1kWh = 3,6. 10° J, 
das Elektronenvolt 
leV 


die Kilokalorie (Internationale  Tafel- 
Kilokalorie vom Jahre 1956), definiert als 


1 keal 


= 1,602. 10-'* J, 


4186,8 J. 


3.3 Die Einheit der Leistung ist das Watt (W), 
IW=1J/s. 
3.4 Die Einheit des Druckes ist das Newton je 
Quadratmeter, 
1 N/m*® = 1 kg m= s~*, 


Als weitere Druckeinheiten kénnen verwen- 
det werden: das Bar, 

1 bar = 105 N/m, 
die Technische Atmospahre (at), 

1 at = 1 kp/em? = 98066,5 N /m?, 


die Physikalische Atmosphare (atm), definiert 
als 
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latm = 101 325 N /m?, 


das Torr, definiert als 


1 
T = —— at 
1 Torr 760 atm 


(= oe Nm? & 133,32 N m?). 
760 

4. Es sind Gréssengleichungen zu bevorzugen. 
Bei der Auswahl der in diese Gleichungen einzu- 
fihrenden Gréssen ist zu beachten : 

4.1 Als Mass der Menge ist ihre Masse, nicht ihr 
Gewicht (Masse Fallbeschleunigung) zu 
benutzen. 

Dies bedeutet insbesondere : 


mal 


4.11 Die spezitischen Gréssen der Warmelehre 
sind auf die Masse, nicht auf das Gewicht zu 
bezichen. 

4.12 Die Wichte (auch spezifisches Gewicht 
genannt) ist zu vermeiden und die Dichte zu 
benutzen. 

4.2 
so zu definieren, dass die Gréssengleichungen in 


Elektrische und magnetische Gréssen sind 


rationaler Schreibweise erscheinen. 

4.3 Energie, Arbeit und Warmemenge sind 
Gréssen gleicher Art. Das Warmeaquivalent tritt 
demnach in Gréssengleichungen nicht mehr auf. 

In Punkt 1 der Empfehlungen wird eine friihere 
Stellungnahme aufgehoben, die der weiteren 
Entwicklung im Wege stand; hierauf sei nicht 
weiter eingegangen. Die tibrigen Punkte lassen 
sich deutlich nach zwei Gesichtspunkten einteilen, 
namlich nach Einheiten (Punkt 2 und 3) und nach 
Gréssen (Punkt 4). Diese werden im folgenden 


getrennt behandelt. 


38. Dik EmMPrOHLENBN EINHEITEN 


Zu Punkt 2 

Die lange Diskussion tiber die Benennung der 
Krafteinheit hat ihr 
vorlaufiges Ende dadurch gefunden, dass das 
Technische Kommitee ISO,/TC 12 zwei Namen als 
gleichberechtigt anerkannte [2] : Kilopond (kp) und 
kilogramme-force (kgf). Eine Einigung auf eine 
einzige Bezeichnung gelang nicht, jedoch sollen 
die beiden Symbole kp und kgf die Vielzahl der 
bisherigen ersetzen, wie kg,,, kg, kg’, kg*, Kg, kG, 


metrischen  technischen 


Die IUPAP schloss sich dem Beschluss von 
ISO/TC 12 an. 

In Deutschland wird das Kilopond, das man dort 
um 1939 vorgeschlagen hatte, im Schulunterricht 
gebraucht, auch verwendet es der Deutsche 
Normenausschuss seit 1955 in seinen Norm- 
blattern. In Osterreich ist das Kilopond gesetz- 
liche Einheit und wird auch in das kiinftige 
deutsche Einheitengesetz aufgenommen werden, 


das man zur Zeit vorbereitet.* 


Zu Punkt 3 


Hier ist der Kernpunkt der Empfehlungen 
enthalten, namlich die Bevorzugung des Inter- 
nationalen Einheitensystems. Die Urheber dieser 
Empfehlung waren sich véllig dariiber klar, dass 
man mit einem einzigen koharenten Einheiten- 
system niemals auskommt, da man damit haufig 


sehr unbequeme Zahlenwerte erhalten wiirde. 
Hier bietet sich der Ausweg an, dezimale Teile 
und Vielfache Einheiten zu 
verwenden, fiir die von 10°" bis 10" die wohlbe- 
kannten Vorsatzsilben mit ihren Abkirzungen 
existieren (Millt fir 10-*, Kilo fiir 10° usw.) 
Manche Einheiten dieser Art haben eigene Namen, 
z.B. die Tonne (t) fir 10® kg. Allerdings sind der 
Grundsatzes 


der koharenten 


konsequentien Anwendung dieses 
dadurch Grenzen gesetzt, dass die Zeiteinheiten 
Minute und Stunde nicht dezimal aus der Grund- 
einheit Sekunde abgeleitet sind. 

Die Empfehlung bedeutet doch aber zugleich 
den Vorschlag, nichtkoharente und nicht dezimal 
abgeleitete Einheiten in Zukunft (mit der genann- 
ten Einschrankung) nach und nach fallen zu 
lassen, sofern nicht sehr triftige Griinde fir 
ihre Weiterverwendung bestehen. Fir den Bereich 
der Thermodynamik handelt es sich dabei 
vorzugsweise um die Druckeinheiten at und atm 
und um die Energieeinheit keal. Die technische 
und die physikalische Atmosphare waren durch 
das Bar zu ersetzen, das bereits als Millibar in der 

* Fir den Bereich der Deutschen Demokratischen 
Republik ist die “ Verordnung iiber die physikalisch- 
technischen Einheiten”’ kiirzlich in Kraft getreten 
(Gesetzblatt Teil I. Nr, 56 vom 6. September 1958, 5S. 
647-649). Hierzu gehdrt die “ Tafel der gesetzlichen 
Kinheiten.”’ Mitteilungsblatt Nr. 149 des Deutchen 
Amtes fiir Mass und Gewicht vom 8, Oktober 1958. 
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Meteorologie verwendet wird. An die Stelle der 
Kilokalorie soll das Kilojoule treten. © Naturlich 
braucht eine solche Umstellung cine verniinftige 
Ubergangszeit, weshalb die EKmpfehlung auch 
einige heute tibliche Einheiten aufzahlt. 

Es ist cine sehr bestechende Eigenschaft des 
Internationalen Einheitensystems, dass man mit 
ihm die verschiedenen Energiearten in der gleichen 
Einheit ausdriicken kann. Das wird aus folgender 
Einheitengleichung deutlich : 


N.m kgm’ st = J Wis (1) 


Die cinzelnen Terme dieser Gleichung kann man, 
obwohl sie per detinitionem alle dasselbe bedeuten. 


als Einheiten folgender Energicarten auffassen : 


Arbeit (Kraft mal Weg), 


kinetische Energie (Masse mal Quadrat der 
Geschwindigkeit ) 
Warmemenge, 


elektrische Enervie. 


Es gibt keinen verniinftigven Grund dafir. 
mechanische Knergn nnurin kpm, Warmemengen 
in keal und elektrische Energien in kWh anzugeben 
und dadurch eine Vielzahl von Umrechnungsfak- 
toren notwendig zu machen. 

Die Anwendung des Internationalen Eimheiten- 
systems bietet auch cine Moglichkeit, zwischen 
metrischen und englisch-amerikanischen Einheiten 
zu vermitteln. Weder die Kilokalorie (keal) noch 
die British thermal unit (B.t.u.) sind in irgend 
Der Ersatz 


neue st 


einem Einheitensystem koharent. 
Einheit 
erfahrungsgemass umso leichter durchzufthren, 


einer gewohnten durch eine 
je weniger die Umrechnungsfaktoren von Eins 


abweichen. Da mehrere internationale Gremien 
bereits verabredet haben, sich im internen Verkehr 
nur des Internationalen Ejinheitensystems zu 
bedienen, ist die folgende Zusammenstellung der 
Umrechnungsfaktoren fiir die Druck und Ener- 


eiceinheiten von Interesse. 


l at 

latm 

1 bf in? 

1 kealyy 

1 (Btu. pp 


0.980665 bar (genau) 
1.01325 bar (genau) 
0,0689476 bar 
4.1868 kJ (genau) 
1.05506 kJ 


Dabei ist Ibf dic Abkurzung fur die Kraftemheit 


pound-foree. Unter kealpp und (B.t.u.)pp sind die 
Tafelkilokalorie inter- 
national table British thermal unit zu verstehen, 
Internationale Dampftafel- 


Internationale und die 


die durch die 5. 
konferenz (London 1956) angenommen sind und 
nach threr Definition nicht mehr auf den Stoffeigen- 
schaften von Wasser basieren. 

Die 


Ubergang von at oder atm auf bar die Zahlenwerte 


Zusammenstellung zeigt, dass beim 
Dasselbe gilt fiir den 
Wesentliche 


Umgewoéhnungen waren beim Ersatz von Ibf in? 


fast ungeindert bleiben. 


Ubergang von (B.t.u.))p auf kJ. 


(auch p.s.i. genannt) durch bar und von (keal)py 
durch kJ Falle 
nebenbei der Vorteil, dass die spezitische Warmeka- 
pazitat von Luft etwa den Zahlenwert Eins erhalt 
(bei OC und 1 bar ist Cy 1,00, kJ kg gerd). In 
mehreren Landern (Schweden, Schweiz) gibt es 


notig. Im letzteren ensteht 


heute schon Wasserdampftafeln, die im den 


Kinheiten bar und Joule abgefasst sind. 


t GROSSEN UND 


(ROSSENGLEICHUNGEN 
Die Schwierigkeiten, die aus der Verwendung 
Einheiten 


werden 


verschiedener und Einheitensysteme 


entstehen, wesentlich gemildert, wenn 
man die Formelzeichen nicht als Zahlen sondern 
als physikalische Grossen auffasst. Eine Gleichung 
zwischen solchen Gréssen, cine  sogenannte 
Grossengleichung, beschreibt cin Naturgesetz in 
seiner allgemeinsten Form ohne das Beiwerk von 
Umrechnungszahlen, Die Gréssengleichung gilt 
unabhangig von der Wahl bestimmter Einheiten, 
sie ist auch gegen den spateren Wechsel von 
Kinheiten invariant und legt den Benutzer nicht 


auf die Verwendung koharenter Einheiten fest 
[3]. [4]. [5], [6). 
Zu Punkt 4 

Da man mit Hilfe bestimmter cinfacher Regeln 
Kinheiten 
umformen kann, ist es zweckmiassig, die Gréssen- 


eme Grodssengleichung fir beliebige 
gleichungen stets als Ausgangsgleichungen an- 
zugeben. Ihre Verwendung wird daher durch den 
Vpl Der Vorteil der 


Groéssengleichungen wirkt sich erst dann voll aus, 


besonders —empfohlen. 


wenn sie fiir den Ingenicur in der gleichen Form 


geschrieben werden wie fiir den  Physiker, 
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Chemiker, Physikochemiker usw. Das setzt aber 
bestimmte gemeinsame Regeln fiir die Definition 
der abgeleiteten Gréssen voraus, auf die daher in 
den letzten Teilen der VDI-Empfehlungen ein- 
gegangen wird. 

Zu Punkt 4.1. Natirlich braucht die Definition 
abgeleiter Gréssen nicht von der Wahl irgend- 
welcher Einheiten 
abzuhingen. Trotzdem kann gelegentlich die Aus- 


oder Grundeinbeiten 
wahl der Grundeinheiten auch auf die Groéssen- 
definition Ein solcher Fall 
offenbar beim metrischen technischen Einheiten- 


einwirken. liegt 
system vor. Weil die Ingenieure die Krafteinheit 
Kilopond (kilogramme-force) zur Grundeinheit 
und damit die Kraft zur Grundgrésse gewahilt 
hatten, bevorzugten sie auch das Gewicht (im 
Sinne einer Kraft) als Mass fiir die Menge und 
zugleich als Bezugsgrésse fiir spezitischen Grossen 
Dabei 
nicht um das értliche Gewicht G, sondern um das 


der Thermodynamik. handelte es sich 

Normgewicht G,. das mit der Masse m durch die 
Gleichung 

G, (2 

der g,, Norm-Fallbe- 


9.80665 m_ s*). 


zusammenhingt, in cdlie 
schleunigung bedeutet (g,, 

Nehmen wir wieder die spezitische Warmeka- 
pazitat von Luft als Beispiel, so erhalten wir 
folgende beiden Gleichungen : 
fir die Ingenicure 
kealp, 
0,240 se 

kp grd 

und fiir die Physiker 
kealpy 


m ky grad 


C 
c 0,240 


p (4) 
wobei C die Warmekapazitit bedeutet. Obwoll 


c,” und ¢,. Gréssen verschiedener Art sind, haben 


v p 


die Zahlenwerte den gleichen Betrag. Zwischen 


den beiden Gréssen besteht die Gréssengleichung 
(5) 

und zwischen kg und kp die Einheitengleichung 
kp = kg. g, (6) 

Ublicherweise haben die Ingenieure bisher die 


Grosse ¢,' mit ¢, bezeichnet und die Krafteinheit 


kp mit kg. Sie schrieben also die Gleichung 


keal py 


0,240 (7) 


kg grd ; 
die dusserlich mit Gl. (4) identisch war, in 
Wirklichkeit aber Gl. (3) bedeuten sollte. Diese 
laxe Schreibweise, die manches Missverstandnis 
hervorrief [6], [7], [8], erleichtert jetzt die durch 
VDI emfohlene Umstellung (Punkt 4.11). 
Es ergibt sich das Kuriosum, dass in der techni- 


den 


schen Literatur bisherige Gleichungen nach Art von 
Gleichung (7) stehen bleiben kénnen, obwohl die 
Grosse auf der linken Seite ihre Detinition geandert 
hat und die Einheit auf der rechten Seite zu einer 
anderen Gréssenart gehort. 

Die im metrischen technischen Einheitensystem 
kohirente Masseneinheit (kp s* m — 9,80665 kg) 
hatte in Deutschland nie einen eigenen Namen. 
Da die spezitischen Grossen, wie oben erwalnt, 
bisher auf das Normgewicht (in kp) bezogen 
wurden, kam diese Masseneinheit als Bezugsein- 
heit Frage. 
Umstand kénnen bei der empfohlenen Umstellung 


nicht in Durch diesen gliicklichen 
die Zahlenwerte der spezifschen Groéssen ungein- 
dert bleiben. 

Auch das spezifische Volumen v wurde bisher 
durch den Ingenieur anders als durch den Physiker 


definiert. Fir den Ingenieur war 


i l/y 


y als Wichte (spezitisches Gewicht), wobei 
Wichte bei Norm-Fallbe- 


war, fir den 


mit 


eigentlich die der 


schleunigung g, gemeint wahrend 


Physiker v die reziproke Dichte p bedeutete : 
v l p 


Der Ingenieur bezeichnete bisher ve’ mit v und 
unterschied nicht zwischen g und g,, wodurch 
manche Gréssengleichung eine fast unverstandliche 
Form erhielt [6], [7], [8], 
keiten zu beseitigen, empfiehlt der VDJ im Punkt 
+. 12, Wichte 
benutzen. Das spezifische Volumen soll nach 
Punkt 4.11 0 
kraftproblemen wirde man die 6rtliche Wichte 


Um diese Mehrdeutig- 


die iiberhaupt nicht mehr zu 


1 p detfiniert werden. In Schwer- 
besser durch p g ersetzen. 

Die praktische Folge dieser Emphfehlungen 
wiirde sein, dass die Ortliche Fallbeschleunigung 


g nur noch in Gleichungen auftritt, die Schwer- 


137 





U. Grieut. 


kraftprobleme beschreiben. Die Norm-Fallbe- 
schleunigung tritt iberhaupt nicht mehr auf, die 
Wichte (spezitisches Gewicht) ist entbehrlich. 

Der Punkt 4.2 Empfehlungen, der sich mit 
elektrischen Grissen beschiftigt, sei hier tiber- 
gangen. 

Zu Punkt 4.3. 
Energie, Arbeit 


Es ist seit langen bekannt, dass 


und Warmemenge Gréssen 


gleicher Gréssenart sind. Das bedeutet aber, dass 


jede zugieich als Einheit der 
Warmemenge aufgefasst werden kann, z.B. J, erg. 
kWh, kp m, keal, B.t.u., eV. In der technischen 
Literatur herrschte vielfach die Ansicht vor, dass 


mechanische Arbeit nur in kp m, Warmemenge 


Energieceinheit 


nur in keal gemessen werden diirfte, weshalb die 
Umrechnungszahl zwischen diesen beiden Einheiten 


kealpy $26,935 kp m, 


mit einem eigenen Symbol (A) in die entsprechen- 


den Gleichungen aufgenommen und mit der 


besonderen Bezeichnung 


‘ 


*mechanisches Warme- 
aquivalent ” belegt wurde. Da aber irgend eine 
Einheitenvorschrift nicht mit dem Charakter der 
Gréssengleichung vertraglich ist, darf in Grossen 
keine Umrechnungszahl 


gleichungen — auch 


zwischen Einheiten erscheinen, was in Punkt 4.3 
der VDI-Empfehlungen ausdriicklich festgestellt 
wird. Auch verdient die Zahl 426,935 sowenig 
einen eigenen Namen oder ein eigenes Symbol, wie 
irgend eine andere Umrechnungszahl zwischen 
Einheiten. 


5. SCHLUSSBEMERKUNG 


Die Empfehlungen des VDI, die vorstehend 
kommentiert wurden, werden seit einigen Monaten 
in einem Teil der deutschen technischen Literatur 
angewendet. Obwohl die Beobachtungszeit noch 
sehr kurz ist, lasst sich schon heute erkennen, 
dass das Interesse der Fachwelt sehr gross ist und 
dass ernsthafte Schwierigkeiten nicht aufgetreten 
Die 


uberwiegen bei weitem die negativen. 


sind, positiven’ Kommentare der Leser 

Die Autoren dieser Empfehlungen hatten sich 
das Zicl Verstandigung 
zwischen Ingenicuren und Naturwissenschaftlern 
Die Chancen, die hierzu in der 


und 


gesetzt, eine bessere 


zu erreichen. 


Benutzung von Gréssengleichungen von 


einheitlichen Einheiten liegen, sollen voll aus- 
genutzt werden. Dass dabei auch alte Gewolhn- 
miissen, ist unver- 


heiten aufgegeben werden 


meidlich. 
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The influence of the initial filtration velocity and of vibrations on the 
resistance of a polystyrene filtercake 


P. J. Bakker, P. M. Heertyes and J. L. Hrsovu 


Laboratory of Chemical Engineering of the Technical University, Delft. 
(Received 8 March 1958 ; in revised form 27, August 1958) 


Abstract—The influence of the initial filtration velocity and of vibrations on the resistance 
of a filtercake of polystyrene spheres of known diameter and formed over three different well- 
defined filters has been studied. 

The influence of the initial filtration velocity could be expressed mathematically, showing 
an increase of a cake resistance zero at an infinite velocity, to a maximum value at a velocity 
zero. The parameters in this equation depend besides on the system on the particle size. 

The influence of vibrations on the cake resistance has been studied qualitatively. In general 
it has been observed that the mean resistance increases by vibration and with increasing amplitude 
and frequency. The influence is larger if the cake has a high resistance and is not very stable. 
For the case of a well-defined vibration, a mathematical expression for the relative cake resistance 
as a function of the frequency could be derived. 

By successive filtration, filtration under vibration and re-tiltration, interesting phenomena 
have been observed. 

A qualitative explanation of the phenomena, based on the structure of a cake as a function 
of the initial filtration velocity and based again on the structure of foregoing layers, has been 
given. 


Résumé —-L ‘influence de la vitesse initiale de filtration et de vibrations sur la résistance d'un 
gateau de sphéres de polystyréne de diamétre connue forme sur trois filtres différents bien définis, 
a été étudiée. 

L’influence de la vitesse initiale de filtration a pu étre exprimée mathématiquement, 
montrant une augmentation de la résistance du gateau, nulle pour une vitesse infinie, & une 
valeur maximum pour une vitesse nulle. Ces paramétres de cette équation dépendent du systéme 
et de la grandeur des particles. 

L/influence des vibrations sur la résistance du gateau a été étudiée qualitativement. En général 
il a été observé que les vibrations augmentent la résistance moyenne et ceci d'autant plus que 
amplitude et la fréquence sont plus grandes, L’influence est plus grande si le gateau a une 
grande résistance et s'il n'est pas trés stable. Pour le cas d'une vibration bien définie, une ex- 
pression mathématique donnant la résistance relativé du gateau en fonction de la fréquence 
a pu étre trouvée. 

Dans le cas d'une filtration suivie de filtration avec vibration puis de nouveau par une 
filtration simple, des observations intéressantes ont pu étre faites. 

Une explication qualitative du phénoméne, fondée sur la structure du gateau en tant que 
fonction de la vitesse initiale de filtration et sur la structure des couches précédents a été donnée. 


Zusammenfassung—Der Einfluss der anfinglichen Filtrationsgeschwindigkeit und von 
Schwingungen auf den Widerstand eines Filterkuchens aus Polystyrolkugeln bekannten 
Durchmessers wurde untersucht, wobei der Filterkuchen iiber drei wohldefinierten Filtern 
gebildet wurde. 

Der Einfluss der anfinglichen Filtrationsgeschwindigkeit konnte mathematisch ausgedriickt 
werden ; der Kuchenwiderstand Null bei unendlicher Geschwindigkeit erreichte einen Maximalwert 
bei einer Geschwindigkeit Null. Die Parameter dieser Gleichung hingen ausserdem von der 
Teilchengrdsse ab. 

Der Einfluss von Schwingungen auf den Kuchenwiderstand wurde qualitativ untersucht. 
Im allgemeinen wurde beobachtet, dass der mittlere Widerstand infolge von Schwingungen und 
mit steigender Amplitude und Frequenz anstieg. Der Einfluss ist grésser, wenn der Kuchen 
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einen hohen Widerstand hat und nicht sehr fest ist. Fiir den Fall einer wohldefinierten Schwingung 


konnte ein mathematischer Ausdruck fiir den relativen Kuchenwiderstand als Funktion der 


Frequenz abgeleitet werden. 


Wenn einer geowOdhnlichen Filtration eine solche mit Schwingung und dann wieder cine 


vewOhnliche Filtration folgten, konnten interessante Beobachtungen gemacht werden. 


Die Erscheinungen konnten qualitativ gedeutet werden, wobei man von der Struktur des 


Kuchens als Funktion der anfiinglichen Filtrationsgeschwindigkeit und von der Struktur der 


vorhergehenden Schichten ausying. 


INTRODUCTION 


Ix filtration, vibrations often have a strong effect 
on the filtration velocity. In former publications 
[1], [2] and [3] precautions have been described 
which had to be taken against an influence of 
vibrations on the results, 

Because of the technical importance of this 
phenomenon it has been undertaken to study the 
the 
To facilitate a possible explanation, 


influence of vibrations on resistance of a 
filtercake. 
a well-defined system has been chosen, viz. 
slurries of spherical particles of polystyvrene* in 
water, filtered over well-detined filters. 

The use of these systems enabled us to study 
more quantitatively than was previously possible 
[2] the influence of the initial filtration velocity 
and (or) of the filter resistance on the packing of 
the filtercake. 

EXPERIMENTAI 


used was almost 


., 
rhe 


another set-up deseribed earlier {1}. Nevertheless 


apparatus equal to 
a sketch (see Fig. 1) and a short description of its 
use will be given. 

The slurry to be filtered is prepared and stored 
in a vessel A (30 1) provided with a stirrer and a 
heating coil controlled by a contact thermometer. 
Via a membrane valve the slurry flows into the 
actual filtration apparatus B. This apparatus 
(see detail) consists of a perforated brass plate, 
on top of which the filtercloth is placed, a brass 
funnel and a cylindrical plexiglass container 
provided with an overflow. 

The filtrate formed is assembled in the con- 
tainer C, the filtrate volume can be measured by 


means of the calibrated burette D. Vacuum is 


*The Werke 
Germany, had the courtesy 


Hliils ** at 
to put 


* Chemische Recklinghausen, 
a large amount of 
polystyrene pearls of wide size spectrum (30-3000 ,) at 
our disposal, Our thanks to this firm for their support 


may be expressed here, 


applied to C and D by means of the water-jet 
pump H. The big vessel F connected with the 
vacuum system serves as a buffer, the needle-valve 
G enables the introduction of a small but constant 
amount of air into the system. In this way a 
constant vacuum (read on manometer E) can be 
maintained under the filter, ensuring a constant 
filtration pressure (vacuum + hydrostatic water- 
head). The suspension overflowing from B into I 
via J and K, by 
After each the 
filtrate is assembled in L, pumped to N by M, 


is periodically returned to <A, 
means of compressed air. run 
pretiltered over a small filter P to remove algae, 
fine dust particles, ete., and stored in Q. This 
prefiltered liquid is used to prepare the slurries in 
A to which it is transported with the aid of com- 
pressed air. As liquid water has been used, small 
amounts of a wetting agent (Teepol or Triton X 
100) have been added to the water to promote 
monodispersivity. The water temperature was 
8-937 10%, ; 1-045). 

With the apparatus as described filtration 


er ’ 

25 ¢ (water polystyrene 
pressures equal to 17-10* (head of ~ 170em of 
water) have been used. 
smaller filtration 


For a much 


(s-108 


pressure 
~s8em of water) the apparatus was 
changed into a gravity filter. The vacuum system 
was disconnected and the filtrate flowed via an 
overtiow into the container C, No details will be 
given. 

filter “cloths” have been used: a 
nvlon filter (nylon : 


nickel 


Three 
$,000 openings em*), and 
two 


filters (Nil: 729 circular openings 
(6 = 250 u) em®; Nill: 2704 circular openings 
(d 75 jp) em?, 
graded in fractions by sieving. The following sieve 
125-150 pw, 150-175 p, 


The polystyrene mixture was 


fractions have been used : 
175-210 M 350) 420 M. 500-600 Mt. 
710--1,000 fe. 1,200.-1,400 pe 

1,700.-2,000 ¢ and 2,000-2,400 ju. 


600-710 ja, 
1,400-1,700 ju, 
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DETAIL B 
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Fig. 1. Sketch of the filtration set-up. 
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The filter was mounted in a frame. Vibrations 
were applied to this frame. Two different devices 
were used. The first consisted ot a wheel rotating 
around a horizontal axis. To this wheel pieces 
of rubber-tubing were attached in which small 
iron rods were fastened. The rubber ends touched 
the bearing bars of the filter support with a 
frequency proportional to the number of revolu- 
tions of the wheel. The amount of impacts 
per minute on the bars will be called the frequency 
of vibration (F’). 

The second device consisted of a constant speed 
motor, coupled to a variator which enabled to 
regulate the speed between 500 and 3,000 rev min 
(F). The variator was coupled to an excentric 
allowing for an amplitude varying from 0-5 mm. 
The vibration created this way was of a sinusoid 
form. 

The magnitudes measured during a run were 
dV* z@and Vasa function of the time of filtration 
@. Measurements were carried out in the laminar 
regime of flow. The specific resistance of the cake 
r, the resistance of the cake R,, the resistance of 
the filter R,, the total resistance of cake and filter 
R and the initial filtration velocity U,. were 
detined and calculated by means of : 

1 dV_ 1 AP 
A » R 
LAP 


| 
= Ky 


A dé Vad 
R= R, + R, 


_ AdR, AAR 


ec dV cdV 


Measurements were carried out at a constant 


filtration pressure. Unless extraordinary effects 
occurred straight R-V lines were obtained, from 
the slope of which r was found. This r is a mean 
resistance. 
RESULTS 
(with Borum, Kroon and Vernnor) 


(a) The influence of the filtration velocity on the 
cake resistance 


It has been described earlier [2] that the cake 


*See nomenclature. 


resistance as far as filtration variables is concerned 
depends on the concentration of the slurry and on 
the initial filtration velocity U,, or if the filtration 
pressure is constant, on the resistance R, of the 
filter. An increase of r with decreasing U, has 
been observed. No mathematical relation could 
be found, because of a too large scattering of the 
results. This was in all probability caused by 
the not too well defined system used. 

The better defined system used here has enabled 
us to find for this system a mathematical relation 
between r and U, (or R,), which takes the form 
of an exponential function : 


r rq CX} (3: J’ 
- r,exp — (22 
. Ry 


r=r,exp — ( 


Combination of equations (2), (5) and (6) gives 
b, b, = VP/». In these equations r, signifies the 
specific resistance of a cake formed at an initial 
filtration velocity zero (or for a finite filtration 
pressure, at an infinite resistance), in which case 
the filtereake will be of the densest packing 
possible. Moreover, the packing of a filtercake 
formed over a filter which resistance approaches 
zero will approach that of the slurry. For the 
dilute suspensions used the resistance of such 
a cake will be very small. It has been assumed 
that the cake resistance for R, = 0, will be zero. 
The experiments carried out do not contradict 
the last supposition. 

In Table 1 the values of ry, b,, b, and d as 
calculated for six series of polystyrene in water 
been collected. The conditions of the 
experiments were: concentration = 6-66 x 10-3, 
AP = 163 x 10, the R, of the nylon cloth used 
has been changed from 20 « 10® to 285 « 104, 

In Fig. 2 the results for the series 4 and 6 are 
represented graphically, the dots are the measured 
data, the crosses are calculated with the aid of 
equation (5). 

In Fig. 3, r has been plotted as a function of 
d, (arithmetic mean) at constant values of R;. 
The one system investigated so far does not yet 
permit us to draw any conclusions as to the 
significance of the magnitudes introduced. It 


have 
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Table 1. 





Series no. | d, range* in yu 


Pp 
108 
108 
108 
108 
10° 
106 


110 
149 
91 
114 » 
135 x 


162 & 


125 
175 
350 
500-600 
600-710 
710-1000 


150 
210 
420 


< 108 
x 108 
x 108 
x 108 
51 x 108 
52 x 106 


0-126 
0-193 
0-610 
0-572 
0-359 


0-352 








*The d, range given is that of the corresponding sieve openings. 


may suflice to note that apparently these magni- 
For the 
polystyrene cakes investigated the lowest values 


tudes depend on the particle size. 


of r, and b, and the highest values of d and 6 
are found for particles of 350-420 ». The change 
with particle size is gradual with the exception 
of r, for the small-sized fractions. It has to be 
observed, however, that the inaccuracy in the 
calculated value of r, may be large. On the other 
hand, the values as found lie within the actual 
possibilities. This conclusion can be drawn from 
percolation data on consolidated polystyrene 
cakes, obtained by introducing polystyrene spheres 
into the filter apparatus containing water, con- 
solidated by careful and 


prolonged tapping 





150 1% 


+— 
> 


++ . + + . + 


| = : 
1710-1000 p }-—+ | 
| | 











Fig. 2. 


against the wall of the filter. This has been 
continued until the cake has the highest possible 
packing density under these conditions. 

Some of the resistances found with these cakes 


were : 


1280 x 108; 
865 108 ; 
620 108; 
413 108, 


1200-1400 ju, 7 
= 1400-1700 py, r 
1700-2000 ju, 7 
2000-2400 ju, 7 


Analogous experiments with glass spheres (see 
later) gave the results that for: 


= 362 x 108; 
206 « 108, 


d,, = 1400-2000 p, 7 


Pp 


d,, = 2000-4000 yu, r 





6 
15010 


— 


- 500-600 p at | 


~~ + ——4—-+- + ne 














010° 


The specific resistance of a cake r as a function of the 


resistance of the filter Ry for two sieve fractions of polystyrene. 
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Fig. 3. The specific resistance of cakes of polystyrene : 
a function of the particle diameter ¢,,. 
showed excellent reproducibility, 


The 


demonstrated by the fact that the total resistance 


results 


was linear with the amount of particles. 

From the above it can be concluded that the 
r, Values as found are possible. The intluence of 
the particle diameter as found for the cake 
resistance of the consolidated cakes has not been 
found for the r, the 
resistance of cakes formed by filtration at the 


values which represent 
highest packing density possible under filtration 
conditions, Of course this does not necessarily 
need to be the case, because the mechanism of the 
cake, built-up for both types of cake, is in fact 
quite different. 

From Fig. 3 it can be concluded that the general 
tendency is that at constant R,. r inereases with 
d,,. although maxima and minima occur. 

The differences observed in the resistances of 
cakes of polystyrene- and glass-spheres seem to 
the the 
surface-tension solid—liquid on the cake-resistance. 
Experiments different 
quantities of Teepol, presented in Fig. 4* as the 


indicate existence of an intluence of 


with water containing 
cake resistances, measured at constant R,. for 
the fraction 175-210» and AP — 16-3 « 104, as 
the tension liquid- air 


(eq) show such an influence if as may reasonably 


a function of surface 


*The number near the points represent the amount of 
Teepol added in ml 301. of water. 








' 
| 
0 





+. The specific resistance of a cake in liquids of 
different surface tension. 


be assumed the change in o;¢ is accompanied by 
a change in og;. 

It seems worthwhile to investigate what actual 
change in og, causes the observed change in r. 


This has not vet been done. 
(4) The effect of vibrations 
(() The filter. 


tions do not cause a change in the filter resistance, 


It has been found that vibra- 


except in those cases where air is trapped under 
the filter, the openings are small (~ 75) and 
BOO0), 
however, the observed increase in resistance is 


the frequency is high (F In these cases, 


small and of the order of some per cent. 
(i) The cake. 
the resistance of cakes. 


Vibrations in general increase 
The effect increases with 
This 


holds for vibrating while filtering or while per- 


increasing amplitude and cake-resistance. 
colating. Reproducible results are not easy to 
obtain especially with particles of large diameter 
and sometimes with cakes of high resistance. 


(ce) Combined filtration and vibration 

In those cases, where vibrations were applied 
from the very beginning of the filtration, linear 
R Vv 
to determine the resistance rp, of the vibrated 
cake. 
expressed in a relative increase of the cake 
(Tp 
of a cake formed without vibration. 


relations were obtained, which enabled 
The effect of the vibrations could then be 


resistance Ar r)/r, if r is the resistance 

It has been found that the relation between 
Ar and F at a constant amplitude again can be 
expressed as an exponential function. Including 
that for F — 0, Ar 
fits the results : 


0, the following equation 


Ar — aexp (*) (7) 
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Table 2. 





d,, range 
(m) 


Amplitude 
Cone. (mm) 'P-+a 
140 
146 
1-50 


1200 
1100 
1000 


106 
106 
106 


17-0 
16-8 
171 


106 


350) 420 


O74 106 


1800 
850 


1-36 
0-61 


106 
106 


1-70 


1-95 


OW 1000 2-92 0-89 106 


S00 GOO O-t4 


900 108 


650 


106 
108 


2-68 
1-81 


1-00 
1-04 





*In this case the amplitude could not be measured, but nevertheless the cake was vibrated. 


The value of a signifies the maximum relative 
increase possible at a certain amplitude, 

In ‘Table 2 the calculated values of a, b and d, 
obtained with experiments where vibrations were 
applied from the very beginning, with a nylon 
filter (R, = 2-69 
of 7-8 
have been collected, The maximum (extrapolated) 
0 (rp... (1 
well as the resistance of the cake formed without 


10°), with a pressure difference 
10° and over an F range of 500-3000 


value of rp, for F a)r| as 


vibration (r) has been included. The vibrations 
were applied either on top or at the bottom of the 
filter. 

The relation between r» of a vibrated cake and 
R, can also be expressed by equation (5). How- 
ever, as compared with a non-vibrated cake the 
magnitude of the constants changes. For a sieve 
fraction 175-210 ju, Fi 540, AP 17-4 104, 
c = 6-66 10°-%, R, region of the Dynal filter 
used 45 10° -+ 230 
have been calculated from the results : 


10°,* the following values 





b, 


10° 
10° 


130 
170 


66 


2-3 


without vibration : 
with vibration : 





*The fraction taken was from a different stock as used 
before. of the 
used normally, a Dynal filter cloth was taken because 
this could be cleaned with diluted hydrochloric acid, 


Moreover instead of one three filters 


From these data it can be calculated that the 
value of ry rat Rp — 40 «© 10% is equal to 1, that 
it rises rapidly to 3-1 at R, — 60 « 10° and then 
decreases slowly to 23 at Ry, = 200 « 10% In 
view of the accuracy of the measurements, the 
difference between the ratio 3-1 and 2-3 is not 
significant. 

The conclusion can be drawn that the increase 
of the resistance by vibration is larger if the cake 
is formed at a relatively low initial velocity and 
therefore has a corresponding relatively high 
resistance, showing an increasing influence by 
decreasing the flow velocity, reaching constancy 
below a certain velocity. 

(d) 


Alternating periods of filtration and vibration 


Interesting phenomena are encountered if a 
filtration is started without vibration (period 1), 
if then vibrations are applied either while con- 
tinuing the filtration or while percolating (period 2) 
and if finally a filtration without vibrations is 
carried out (period 3). The phenomena observed 
depend strongly on the type of cake formed in 
period 1, Two limiting cases will be considered. 
If the cake is formed at a high initial velocity 
and therefore has a low resistance, vibration in 
period 2, while filtrating, immediately shows a 
change in the R-V relation. The straight R-V 
line obtained has a larger slope than in the period 
1. Continuing the filtration without vibrations in 
period 3 results again in a straight R- V line, with 
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the same or a smaller slope as measured in period 
1. The cakes thus formed are of a relatively low 
resistance. 

For a cake formed at a low initial filtration 
velocity, therefore with a rather high resistance, 
the behaviour is quite different. As soon as the 
vibration has been started a sharp and, compared 
with the first case, relatively larger increase in R 
occurs. A second difference is that the R-V 
line for the period 2 is not straight, but curved 
downwards. On stopping the vibrations but 
continuing the filtration (period 3) the behaviour 
again shows a difference with the former case. 

First a sharp drop in the total resistance occurs; 
this can be followed by an increase of R. In some 
cases several repetitions take place with a gradual 


The influence of vibrations on the form of the RV line. 


decrease in the differences between the consecutive 
maxima and minima. Finally again a linear 
R-V relation is found if the filtration is prolonged 
long enough. The slope of this line is larger than 
that of the original cake, the limiting case being 
that by extrapolation the line passes the starting 
point of period 1(V Oo, R R,). 

For inbetween cases the phenomena are in- 
between. The smaller the original cake resistance, 
the smaller the drop of R at the end of period 2, 
the more the undulating part of the R-V line 
disappears, the smaller the increase in slope of 
the straight RV line in period 3, as compared 
with that in period 1, It is not necessarily so that 
the drop in R occurs immediately after the 
vibrations have been stopped. 
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To substantiate the above three examples have 
been presented in a graphical form in Fig. 5. 
The carried out for the 


fraction 


filtrations have been 
75-2104, in a 
6-66 10°%, filtered over 
AP = 163 = 10%. The 
period 2 of vibration (F! 
cloth resistance several layers of nylon have 
been used or the cloth was not cleaned thoroughly. 
Line 1 represents the results of a cake with a 


concentration of 
cloth 
arrows mark the 


nylon with 
two 


540). To increase the 


relatively low resistance (note the value of Ry, to 
be found by linear extrapolation of the R-V line 
of period 1). 

Line II is for a cake with a higher resistance. 
This is an inbetween case as described. 
Immediately after vibration is stopped, a normal 
‘ake filtration follows, the cake showing a some- 
what higher resistance than that of the cake formed 
in period 1. 

Line III is for a still higher original cake- 
resistance showing a still sharper increase in R 
on vibration in period 2, followed by the period 
of unrest mentioned. 

The straight R—V line of period 3 has not yet 
appeared here, because the filtration has not 
been prolonged long enough. 

The phenomena just described can be altered 
by uncontrollable outside influences or even by 
small experimental changes. It is for this reason 
that can 


vet be given. 


no rigorous quantitative treatment 


Some of these influences need 
mentioning. 

It can happen that while carrying out a 
filtration the resistance suddenly increases above 
the normal, especially if slurries with large 


particles are used. In most cases these changes 


Table 


could be attributed to some outside vibration, 
such as a resonating water tap. 

In the second place, especially for thin cakes 
and with small driving forces the constancy of 
the slurry inlet is very important. The regulation 
is not easy for small quantities. A too large 
inlet stream causes a whirling-up of the cake 
layers already formed, causing a 
Moreover, with these thin cakes the 


decrease in 
resistance. 
disturbance can cause deblocking as well as block- 
ing upon vibration. It has often been observed 
that vibration of a thin cake was accompanied 
by a decrease in the total resistance, sometimes 
without an increase afterwards. This happened 
more frequently with nickel filters, than with the 
nylon filter. A change of r from 7-2 « 10® to 
1-7 = 10% can occur in such cases. 

This influence has been avoided as much as 
possible, but could not always be prevented. 

In the third place it may be remarked that air 
under the filter causes a small increase in the 
the filter. For nylon cloth R, 
changes from 2-69 « 10, to 2-87 & 10,. 

In view of these facts, only a few data will be 
given. These data (see Table 3) will be presented 
as mean values of a number of experiments of the 
cake resistances before and after the period of 
3000), therefore of period 1 and 


resistence of 


vibration (F 
of period 3. 

In series 1 and 2, carried out with thin cakes, 
always a permanent decrease of R on vibration 
has been observed. In the other series on vibra- 
tion after a small of R an 
followed, to higher value than that of the final 
value at the end of period 1. 


decrease increase 





Number of 


experiments Cone. 


Filter 


| 666 » 


| 
| Ni | WW 
6.66 


nvlon 





10 x 10-4 
10-3 
333 x 10-3 
10-3 
10-3 


r after 


dy range r before 
vibration 


(») vibration 
0-14 « 108 
O14 « 106 
2-78 « 108 
285 « 106 
228 « 105 


0-19 »« 108 
O15 x 108 
2-63 « 108 
2-67 « 108 
1-67 = 108 


- 


— 
& & 


—-—— 
~~ -) «3 
oe Ge Ge 














*During vibration in this series no filtration took place, the cake was percolated. 
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The results given substantiate some of the 
statements made above. 


Discussion 


The key to an understanding of many of the 
results obtained may be the results published 
earlier [2] and repeated here in a more quantita- 
tive manner, on the influence of the initial 
filtration velocity on the building up of the cake. 
For the filters used the very first layers of the 
cake formed over the filter govern the building 
up of the rest of the cake, of course under the 
thick and 
phenomena such as packing compressibility do 


condition that the cake is not very 
not occur. The packing of the very first layers of 
a cake is now governed for a certain filter by the 
initial filtration velocity. At a high velocity, a 
cake is built up of a lowt density and of a low 
resistance ; at a low initial filtration velocity a 
cake of higher density and of greater resistance 
is formed. 

For the less dense-packed cakes, as a result of 
the great flow velocity, large stabilizing forces or 
pressures must exist. The stabilizing forces are 
the 
maximum near the filter and zero at the top of the 
cake. 
density of the packing and therefore the resistance 
of the cake. 

If the balance of forces is disturbed, e.g. by 


counteracted by cake pressure, which is 


The cake pressure tends to increase the 


vibration, an increase of density is to be expected, 
the more so, if the stabilizing forces are small, 
the frequency is high and the amplitude large. 
Let us now, as a limiting example, consider 
the behaviour of a cake of low density and rela- 
tively great stability under flow and vibration. 
Vibrations of low frequency will not have much 
effect on such a cake, however increasing the 
frequency will finally result in a cake consolida- 
tion starting near the filter. If during the vibra- 


tion filtration is continued, the layers on top 


(which are in their original packing) will determine 
the structure of the newly formed layers giving 
these the original top packing density, whilst the 


tThe indication high and low are used here os before 
ina relative sense. All the cakes considered are far removed 
from the densest packing possible. 


consolidated layers near the filter will grow. The 
effect resembles the retarded packing compressi- 
bility as described by Rrerema [5}. 

Therefore a constant, but higher, mean resistance 
will be found for the vibrated cake. If the vibra- 
tion is ended, the top layers will still have the 
corresponding original packing and therefore the 
same or a smaller cake resistance as that of the 
original cake will be found on continued filtration 
for the newly formed cake. 

As a second limiting case, the behaviour of a 
cake formed at low initial filtration velocity of 
high density and low stability will be considered. 
At a lower frequency than in the former case, 
vibrations will affect the cake density, further- 
more the layers affected will be thicker. If it is 
assumed that the whole layer already formed 
will be affected, the results will be twofold. 
First of all in the beginning a rapid increase in 
resistance, which is much larger than in the first 
case, will occur. Because the driving force will 
decrease with time, the relative effect will decrease. 
Therefore the R—V line obtained will be concave to 
the Vo axis. In the second place the top layers of 
the cake, when vibration is continued long enough, 
will be of a higher packing density than that of the 
original cake. A cake formed on these layers by 
normal filtration will therefore have a_ higher 
density than the top layers of the original cake. 
The resulting mean resistance of the cake cannot 
be higher than when the total cake has been 
formed under vibration from the very beginning. 
Therefore the extrapolated R-V line for the last 
filtration in this limiting case must pass through 
the point Y= 0, R= R,. The just 
discussed have been realized, as can be seen from 


cases 


the results presented. 
The the 
phenomena mentioned follows from the above. 


qualitative explanation of other 
When a cake is formed directly under vibration, 
a straight RV line might be expected, and has 
been found. 

The irregular behaviour of vibrated cakes of 
high density and resistance, when vibration is 
stopped, cannot, however, be deduced from the 
above. It is possible that an explanation lies in 
the observation that vibration will cause a state 


of unrest in the cake. It is conceivable that if 
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the vibrations have not yet resulted in a packing 
stable with respect to conditions of flow, stabiliz- 
ing forces and vibration, the particles will adjust 
themselves again to the new conditions of flow. 
This phenomenon wil! be more marked with less 
cakes. It that an 
analogous behaviour have been described earlier 
[4] for cakes consolidated by flow only. 


stable may be remarked 


CONCLUSION 
For the type of relation between the cake 
resistance and the initial filtration velocity an 
found, The 
parameters in the equations derived depend on 


exponential function has been 
the size of the particles and on the concentration. 
The same type of correlation has been found for 
the relation between Ar, the fractional increase 
of the cake the 
frequency of the vibration. Here, for one particle 


resistance by vibration, and 
size besides depending on the initial filtration 
velocity and on the concentration, the parameters 
depend on the amplitude of the vibration, 

the cake 


resistance. The increase in general is larger with 


Vibrations cause an increase of 
original 
resistance. A the 
behaviour of cakes under vibration, and after 


increasing frequency, amplitude and 


possible explanation of 

vibration, has been given, based on the distur- 

bance of the balance between stabilizing forces 
= 

and cake pressure and on the phenomenon that 


the packing of a layer of particles is determined 


by the packing of the foregoing layers, whereas 


for the type of filters used the packing of the 
first layers is determined by the velocity with 
which these particles reach the filter. On these 
principles the results obtained could be explained 
qualitatively. 


NOTATION 
maximum relative increase of resistance 
» 
em 
1 
1 
1 


surface of filter 

parameter rev, min 
parameter em sec™ 
parameter em 
3 


concentration yom 


parameter no dimension 
p (= 10-4em) 
1 


1 


diameter of particle 
frequency rev min 
impacts min 


” 
1 gec-2 


impacts pro minute 

pressure difference over filter gem 

mean specific resistance of a cake of 

one g per cm? filter surface 

cake resistance for Ry x 

cake resistance of a vibrated cake 

relative increase of cake 

resistance no dimension 

resistance of filtercloth plus cake of 

one cm? filter surface 

resistance of filtercloth of one em? filter 

surface em! 
resistance of the cake on em? filter surface em! 

= initial filtration velocity em sec~! 

em? 
1 


volume of filtrate 
surface tension dyne cm 
density 

viscosity 


time 
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Abstract—The factors affecting the flow of granular solids through orifices of different diameter 
and shapes under gravity conditions has been investigated at various heads, Different granular 
materials such as sand, sugar, rape seed, wheat and rice were used and the best equation was 


found to be : 
Ww D, 0-185 
—_— 0-2: (7 
pp Avy (22 D,) Ad 


where W is the weight discharged per second, pp, the bulk density of the solids, A is the flow area 
of the orifice of hydraulic diameter D,d is the average arithmetic screen diameter and A is 
the average shape factor of the particles, 

The effect of head and container diameter was found to be insignificant in the 347 runs used 
in the statistical analysis. 

The equation allows the prediction of flow rates through various orifices with an over-all 
accuracy of + 10 per cent. 


Résumé—Les auteurs ont étudié les facteurs qui affectent écoulement de grains solides, par 
gravitée a différentes hauteurs a travers des orifices de formes et de diamétres différents. Divers 
matériaux granulaires ont été utilisés, tels que sable, sucre, graines de colza, blé et riz et Péquation 
qui convient le mieux est : 

" 0,185 


pp Ay (22 D,) 


D, 
dda, 


= 0,236 | 


dans laquelle W — débit massique par seconde, pp = densité globale des solides, A — laire 


écoulement de Voritice de diamétre hydraulique D,. d,, le diamétre moyen arithmeétique 


du tamis et A = le facteur de forme moyen des particules, 

L’effet de la hauteur et du diamétre du récipient a été jugé insignifiant au cours des 347 
essais d'analyse statistique. 

L’équation permet de prévoir les vitesses d’écoulement a travers divers orifices, avec une 
précision moyenne de + 10 per cent. 


Zusammenfassung Die Kinflussgréssen auf das Fliessen korniger Feststoffe durch Offnungen 
verschiedener Durchmesser und Form unter dem Einfluss der Schwerkraft wurden fiir verinderte 
Druckhéhen untersucht. Verwendet wurden verschiedene kérnige Stoffe wie Sand, Zucker, 
Riibsamen, Weizen und Reis. Als beste Gleichung ergab sich folgende : 


Ww - 2 ( D, _ 
ppAvi(zeD,) A dap 

Hierin ist W der Gewichtsstrom, py die Schiittdichte der Feststoffe, A die durchstroémte 
Flache der Offnung mit dem hydraulischen Durchmesser D,, d,, der arithmetische Mittelwert 
des Siebdurchmessers und A der mittlere Formfaktor der Teilchen. 

Der Einfluss der Druckhéhe und des Behalterdurchmessers ergab sich als unbedeutend in 
den 347 Versuchen, die der statistischen Auswertung zugrunde lagen. 

Die Gleichung erlaubt die Vorausberechnung von Fliessgeschwindigkeiten durch unterschid- 
liche Offnungen mit einer gesamten Genauigkeit von + 10 per cent. 
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INTRODUCTION 
Tnuere has been very little systematic work 
published in the literature in the last 50 years 
upon the flow characteristics of granular solids 
flowing through orifices from hoppers or bins, 
Due to the lack of any formula for 
predicting the flow rate of granular solids the 


reliable 


design of bins or hoppers has been to some extent 
of an empirical nature, It is the purpose of this 
paper to investigate the flow of granular material 
through various orifices from vertical-sided hoppers 
to find the effect of such variables as bulk density, 
head, orifice diameter and shape, particle diameter 
and shape and container diameter upon the 
discharge rate. 

In the past most of the early work upon flow 
of solids has been a by-product of research into 
pressure distributions in hoppers or bins for 
design purposes [1, 2, 3). The first investigation 
into the flow of solids appeared in 1911 by 
Kercnum [4]. He used wheat and showed that 
the flow was independent of the head and varied 
as the cube of the orifice diameter. Kercnum, 
using glass-sided bins, showed that the flow of 
solids was an intermittent flow caused by the 
rapid formation and collapse of domes or bridges 
which built up in the mass of the material. 

Deine and Menrine [5], in 1929, examined the 
flow of granular solids from conical bins or con- 
tainers and proposed an equation of the form : 


mn [34-6 


toe 100 D*? pg _— 
(67-4 + 444 sin + ¢)] {(d ‘D) + 0-130 — 0-161 | 


equation and also does not show how he derived 
it. 

Newron, Dunnam and Simpson [7], in 1945, 
published some results on the flow of catalyst 
pellets from a hopper for feed to a fluidized 
cracking plant. They found for catalyst pellets 
0-1-0-2 in. in diameter that the maximum flow 
rate given by an equation of the form given below : 

W — 8-50 D?* H®®* |b /min (3) 
where H is in ft and D in in, 


They found that the equation (3) was only valid 
if the orifice diameter exceeded six times the pellet 
diameter, as for smaller orifice sizes the flow 
tended to become irregular and was liable to stop. 
They also investigated the flow patterns and 
showed that for heads less than 1} times the 
container diameter that a stationary cone of 
material remained at the walls subtending an angle 
approximately twice the angle of repose of the 
material, 

Brown and Hawks.tey [8], in 1958, made a 
thorough examination of the mechanism of flow 
of solids through orifices and showed that there 
were at least five distinct regions of movement 
near the orifice and that the upper surface of the 
solids did not move until a considerable dilation 
of the packing had occurred. They did not 
however give any equation to explain the flow. 

Rupp [2] and Jenike [3], in 1954, published 


g/min (1) 


with D and d in mm. 


This equation (1) should be valid for flat 
bottomed containers if the value of dé = 180° — 20 
where @ = tan! yw is used as the conical angle. 

Takanasut [6], in 1933, carried out a large 
number of runs using different materials and 
published an equation which he claimed to 
explain the flow of granular solids. 


— k gD 
~ f (p) +a(d/D) 
It will be noted that this equation does not 


involve a head term. TaAKAnasut does not appear 
to have claimed any specific accuracy for his 


V (2) 


papers trying to explain the flow of solids on the 
basis of the pressures set up in a mass of material 
and while they gave interesting theories their 
equations prove to be complex and the factors 
needed difficult to obtain readily. 
FRANKLIN and JOHNSON [9], in 1955, published 
a paper in which they give a proposed flow 
equation of the form : 
W = 
p, D?* 


(6-288 x + 23-16) (d + 1-889) — 44-90 Ib/min (4) 


The equation (4) was derived for circular 





Kh. T. Power 
orifices only under the conditions of particle size 
0-038 to 0-20 in., 7-3 to 676 lb/lb* 
orifice diameters 0-236 to 2-28 in. and an over-all 


density and 


aecuracy of | per cent was claimed. 
They also suggested that the equation for tlow 
rate could be influenced by the following factors, 


(a) Orifice diameter to particle diameter less than 5. 

(b) Column diameter to orifice diameter less than 
6-3. 

(c) Bed height less than one column diameter. 


It will be noted that in equations (1), (2) and 


(4) the factor of coeflicent of friction appears, 


in the case of TAKANASHT’s equation the method of 


calculation of f (4) is obscure and in the 
equations (1) and (4) the determination of the 
coeflicient of kinetic angle of repose ts not easy 
and subject to considerable error, Franxkuin and 
that the 


a characterization 


JoOuUNSON admit coellicient of friction 


seems to be factor for other 


physical properties such as particle shape, 
roughness and effective size and also void fraction 
of the Since it 
been shown by Fowner and Cnopztesxer [10] 
that the 


particle shape, 


bed of packed material. has 


coetlicient of friction is a 


diameter and roughness it was 
deleted from this investigation in favour of the 
factors. It 


that the possible factors that can influence the 


more fundamental was considered 


flow of dry granular solids are as follows : 


(i) The orifice diameter or the equivalent 
hydraulic diameter for orifice shapes 

other than circular 

(ii) The container diameter 

(iii) 

(iv) 

(v) The true density of the material 

(vi) The bulk density of the 
combination of (v), and the void fraction 
of the bed 


The orifice shape 


The true spherical particle diameter 
The head of material above the orifice 


bed which is a 


(vii) 
It will be noted that the fundamental factors 
affecting angle of repose are included in the above 
list either by themselves or in combination with 
others. 

The application of dimensionless analysis to 
the above factors yields an equation of the form : 


and J, 


case of 


function of 
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a? F(R): r) (2): oh (7 ‘ 


which can be written in a more convenient form as 


D, 
a 1J% D f (2) ho) S-(7 *) (6) 


Equation (6) resembles the basic equations 
derived by Deine and Meurine in 1929, except 
that they 
weight of solid which is a reciprocal of equation 


(5). 


It was decided to investigate the equation (6) 


used time for discharge of a given 


by the use of statistical analysis and a series of 


experiments were carried out involving all the 


possible factors. 


EXPERIMENTAL APPARATUS AND 


PRoceDUn! 


Three Pyrex glass cylinders 48 in. in length 


and of inside diameters of 6), 4 and 2] in. were 
selected and clamped vertically in a steel fram- 


work. The 


of 25, 


tubes were marked to indicate heads 
50, 75 and 100em from the base plates. 
The base plate consisted of 18 s.w.g. brass plate 
and was so constructed that an orifice plate could 
be slid into slots in the base plate together with a 
sliding shutter. The orifice plates were made in 


I8s.w.g. brass plate in which the appropriate 


orifices had been cut. The circular orifices were 
cut to size in a lathe but the other shapes such 
were cut and filed 
Due to hand cutting slight deviations 
but the 
areas and perimeters were accurately measured by 
Details of the 


as triangle, square, ellipse ete. 
by hand, 
from the nominal size were unavoidable 
triangulation using a cathetometer, 
orifices are given in Table 1 

The 


rape seed, rice, 


materials used in the tests were sand, 


wheat and sugar. The sand was 
screened into four size fractions, namely — 20 
+- BO, 30 + 40, — 40 + 50 and 50 Tyler 
mesh by the use of a double deck vibrating screen, 
The true specific densities of the materials were 


determined by pyknometers. The shape factor 


for each material was determined by permeability 
measurements on a Lea and Nurse [11] apparatus 
and the shape factor calculated from the Kozeny- 
Carman equation. 


The angle of repose (static) 
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Table 1. Details of orifices used in flow experiments 





Nominal Hydraulic 
Shape of orifice area Actual area Perimeter diameter 
(em?) (cm?) (em) (em) 


Circle 2 2O-157 15 5-066 
Square 20-547 1S-O4: P51 
Rectangle 19-821 18-88 4-197 
Triangle 20-327 3-912 
530 
269 
131 
‘62 
“B57 
339 
“BAT 
O38 
“O81 
“248 
O14 
‘TH2 
73 
“O38 


Circle 
Square 


to te te 


Rectangle 
Triangle 


co — 


Hexagon 


* 
~ 


| Pentagon 
| Ellipse (2:1) 
(3:1) 
(4:1) 
Rectangle (1 :1-4) 
(1:3) 
(1:41) 
(1:47) 
(Tid d)) 


“SOS 
O79 
“D606 
+40 


Circle 
Square 
Rectangle 
Triangle 


“2383 
“135 
O06 O55 
ATG O-006 


Circle 
Square 
Rectangle 
Triangle 





Table 2. Details of physical properties of the materials used in the flow experiments 





Average , True average 
.verage i 
spherical Angle of 


particle repose 


True density screen 
, é shape 
Material Ps particle 
3 factor 
(y mil) size \ 
d (em) Ad d, 


diameter () 


Sand 
20 : 7 O-O661 


30 + 0-0508 


O-o92 

O-O668 
O-O496 
00-0406 
O-4128 
O-3156 
1924 
O- 1095 


= 
te te 


4 4 WO8T2 


= 
~ 
—-t wt oF 


# ¢ 
ww t 


oO 0-0272 
Wheat $ O-3208 
Rice O-2805 
Rape seed 0-1924 
Sugar O-0712 
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for each material also determined. The 


details of the particle properties is given in Table 2. 


was 


A series of experiments were planned to cover 
all possible combinations of orifice size and shape, 
particle size and head and these resulted in a 
total of 470 experiments, 


EXPERIMENTAL RESULTS AND 


CONCLUSIONS 


Statistical correlation of 347 results to determine 
the values of the exponents and constant in 


equation (7) 


Ww .\°/D,\° (DA ay . 
p, AY 2g D, a(’ (2°) (7) (i) (7) 


PR 


revealed that the significance of the terms 


(**) (7) and | 7) was very low and consequently 


could be ignored. 

The value of the oritice shape factor (8) could 
be explained by using the hydraulic diameter of 
the orifice detined as D, — 4 » free area perimeter. 
The final equation evolved based on 247 runs was 


as follows : 


- ) omraS 
<< 4 0-236 | f ,) (8) 
Pp Ay/2g D, d, 

The individual values of the constant in the 
equation (8) for each orifice shape was further 
calculated and is given in Table 3. 

Table 3. Values of the constant in equation (8) for 
various orifice shapes 





No. of runs 


Constant a 


Orifice shape 
Circle ow 
Ellipse 13 
Rectangle 38 
Square 
Pentagon 
Hexagon 
Triangle 


4 
- 


wwe 
hen 
-— 3 ' 


y ty 
= 


ae 


| 
Weighted — 0-236 
Average 


Total 





It is not possible in this paper to give the full 
set of 470 runs performed but the results for a 


random sample of runs is given in Table 4 together 


with the observed and caleulated discharge 
rates. 

The over-all accuracy of equation (8) is + 10 
per cent which is quite good considering the wide 
range of orifice shapes, diameters and particle 
sizes used, 

The results of Table 4 are depicted graphically 
in Fig. 1 where the observed velocity Wp, A 
em. see is plotted against the calculated velocity 
O1SS 

em see, 


of 0-236 y/22 D, | z) 
as 


wD Ss DW 3 4 : 


tty 36 (290, m Sec 


Fie. 1. 
velocity of flow against the calculated velocity for the 
results of Table 4. 


Graph showing the effect of plotting the observed 


Further examination of the complete set of 
470 results reveals that while equation (8) can 
be regarded as sufliciently accurate for most 
design purposes, it is felt that more complex 
expressions involving a parameter expressing 
orifice shape would probably reduce the error to 
even smaller limits. However, the accuracy of 
determination of such parameters as particle 
shape factor, bulk density, angle of repose, ete. is 
of the order of + 10 per cent and the incorpora- 
tion of these factors in any equation will limit 
the over-all accuracy. 

It will be noted that in Table 4 that several 
results are marked with an asterisk. These results 
were carried out at a low ratio of orifice hydraulic 
diameter to mean spherical particle diameter 


154 





The flow of granular solids through orifices 


Table 4. Random set of results used as a basis for equation (7) 





d | Hi PR W obs | 
(Ib/sec) | 


W eale 
(Ib /sec) 


Orifice shape D, 


Materi. ad 
and number (em) faterial 


Run No. toms 


(ib ft?) 


00-0918 


3-22 
0-50 
0-22 
0-080 
3-49 
3-32 
2-88 
0-51 
0-20 
0-080 
0-073 
1-18 
0-189 
0-176 
0-192 
0-197 
0-089 
0-084 
1-14 
0-220 
0-202 
0-205 
1-80 
0-269 
0-280 
0-040 
1-27 
0-254 
0-394 
0-500 


2-94 
9 o” 0-52 
o 88-9 0-24 
90-6 0-079 
SS-5 2-82 
- 3-00 
a 2-82 
0-52 
0-23 
0-083 
0-079 
1-25 
O-185 
0-188 
0-184 
0-197 
0-085 
0-074 
1-13 
0-163 
0-146 
0-154 
1-61 
0-264 
0-270 
0-040 
1-33 
0-207 
0-377 
0-492 


Sand 91-6 
14 | ° : 2-530 
28 . 1-898 ae 
3Y °° | 1-283 % 
07 9° 5-466 00-0496 

103 Square 4511 

105 Triangle 3-912 

116 Square 2-269 

132 Triangle 22 1-400 S 
137 Square 24 1-135 is 
144 Triangle 26 0-996 As 
197 Rectangle 3 4-197 O-1924 
205 Rectangle 2-131 Pa : 
207 Triangle 1-962 
209 Hexagon 2-357 
211 Pentagon 2-339 
216 Square 1-679 | a ‘is 
217 1-566 _ ™ 
233 4197 | Wheat 04128 

*237 2-530 

246 Hexagon 2-357 

*248 Pentagon 2-339 - - 

251 Square 511 0-10905 

261 Rectangle 2-131 

267 Pentagon 10 2-339 

284 Triangle 26 0-996 a = 

2902 Triangle 4 3-912 Rice 0-3156 
*293 Circle 5 2-530 

$24 Rectangle 16 1-762 

4333 Ellipse 13 1.981 


—— | 
2 | Circle 5-066 | 
i 


- 





~~ 
uo oo 








Rectangle 
Rectangle 
Circle 


Sugar 


Sand 0-0406 











D ; , 
*These runs were performed at a low (2*) ratio approaching 6: 


1 and the flow consequently tended to be unstable. 


orifice diameter and shape, that the effect of reducing 
the particle diameter for a particular material 
such as sand was an increase in the flow rate 
to a maximum at a particle size of 65 mesh Tyler 


approaching 6:1. As found by previous workers 
at this value and lower, the flow of granules tends 
to become irregular and sporadic and consequently 
such results have a high error in prediction by 
equation (8). The influence of head upon the flow 
rate has been found to be insignificant in the 
analysis of equation (7) which confirms previous 
work except for Newron, Dunnam and Simpson 
[8]. However, a very slight effect on the flow rate 


approximately and then a decrease with particle 
size. This effect was quite marked for finer sizes 
than 65 mesh as the material lost its “* flowability” 
and tended to flow in lumps or aggregates with 
sudden jerks. Bridging and scaffolding of the 


is noticable at increasing heads but it is debatable 
as to whether this is due to an increase in bulk 
density due to packing of the particles or to a 
head effect in itself. 

It was noticed in the results at constant head, 


orifice became very noticable and with — 100 mesh 
sand the material in the 6in. diameter tube arched 
even with a 4in. diameter circular orifice and 
would only flow when poked. It is apparent 
therefore for steady flow conditions that the 
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particle size has a limiting value for each material 


used. 

The value of the term (p, py) in equation (7) 
proved to be insignificant probably due to the 
small variation of this ratio in the experiments, 
In most cases of packed beds of materials at 
random the bed porosity lies between 30 per cent 
and 50 per cent and the ratio of (p, py) is 1°30 to 
1-50. 
lay between these two extremes and deliberate 


In the experiments carried out the value 


packing of the columns did not give much varia- 
tion in this ratio. It was considered however 
that the incorporation of the bulk density term 
was desirable instead of using the true density 
as it is noticable that any variation of bulk density 
does affect the flow rate, Consequently the value 
of py, was substituted for p, in equation (3). 

The effect of container diameter D_ proved to 
be insignificant upon the flow rate in the experi- 
Generally in the 

d, is never realized and ts 


ments performed. practice 
limiting value of D. 
clearly fantastic. In the experiments the minimum 
value of (D, d,) was 14-1 with the 2] in. tube and 
wheat but no change in the flow rate was detected 


even when the maximum value of 37-1 was used. 


Fow Ler and JJ. 
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NOTATION 

free flow area of orifice 
diameter of circular orifice 
hydraulic or perimetral diameter 

‘ 1 perimeter of oritice 
contamer diameter 
average screen size of particles 
spherical diameter of particles 
cocllicient of 
tan? 


a function involving the 
kinetic angle of friction 0 where yu 
gravitational constant 

head of packing above orifice 
shape factor of particles d, i 
true density of solid 

bulk density of packing 

orifice shape function 

velocity of discharge from orifice 
subtended angle of conical portion of 
contamer 

weight discharged unit time 


a,b. ©, constants 
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Abstract— Experiments have been carried out to determine the variables affecting the angle 
of friction of dry granular materials. It has been found that the angle of friction (@) of solids 
sliding on surfaces depends on the shape factor (/,), specific gravity (s,), and diameter (D) of 
the sliding material and on the roughness («) of the surface. A dimensionless equation has been 
developed statistically and is of the form 


eS 
a € 
3 } ’ Pail 
pe tan ¢ Ri -b (3) C.8, + d 


where a, b, ¢ and d are constants and have values a = 0-2110, b = 03436, ¢ = — 0-O0171 and 
d — 01834. The equation is useful in predicting the sliding angle of dry materials within the 
ranges f, = 0-67 to 1-00, 8, 1-108 to 11-340, D = 0-0272 em to 0-4026 em and ¢ = 0-0272 em 
to O-3134em. The equation developed explains the trends observed by other investigators and 
may be of use in estimating the angle of slope of chutes, slides, ete., and also the angle of repose 
of beds of granules. 


Résumé —Les auteurs ont fait des expériences afin de déterminer les variables qui affectent 
langle de friction de matériaux granulaires secs. Ils ont trouvé que langle de friction (@) de 
solides glissant sur des surfaces, dépend du facteur de forme (f,), du poids spécifique (s,) et du 
diamétre (D) du matériau glissant, et de la rugositeé («) de la surface. Une équation sans dimension 
a été développée statistiquement ; elle est de la forme : 
0,5 
Mm tan 0 = © + (5) +c.8, +d 
Se D 
ou a, b, ¢ et d sont constants et ont pour valeurs : 
a 0.2110; be O34386: cm —OO17T1 et d=—O0,1854 


L/équation est utile pour prévoir langle de glissement de matériaux sec dans les intervalles 
Se = O67 & 1,00, 8 1,103 a 11,340, D 0.0272 em a 0.4026 om et € = 0,0272 em a 0.5134 em. 
Lequation développér explique les tendances observées par d'autres chercheurs, et peut étre 
utilisée pour estimer la pente de chutes, de glissoirs, etc. et aussi Tangle darrét de lits de matériaux 
yranuleés. 


Zusammenfassung Zur Ermittlung der Einflussgréssen auf den Reibungswinkel von 
trockenen, kérnigen Stoffen wurden Versuche angestellt. Danach hing der Reibungswinkel (@) 
der auf Flichen gleitenden Feststoffen ab vom Formfaktor (f,), der relativen Dichte (8), dem 
Durchmesser (D) des gleitenden Materials und der Obertlichenrauhigkeit (€). Kine dimensionslose 
Gleichung wurde auf statistischem Wege abgeleitet, die folgende Form hatte : 
f 0,5 
w= tang =“ +0(<) +c.8, +d 
I D, 

Hierin bedeuten a, b, ¢ und d Konstanten mit den Werten a = 0.2110, b = 0,3436, 
c 0,0171 und d = 0,1834. Die Gleichung erlaubt es, Reibungswinkel trockener Stoffe 
in folgenden Bereichen zu berechnen: f, = 0,67 bis 1,00, 8; = 1,103 bis 11,340, D = 0.0272 em 
bis 0.4026 em und ¢ = 0,0272 cm bis 0,3134em. Die entwickelte Gleichung gibt die Ergebnisse 
anderer Beobachter richtig wieder und kann zur Berechnung der Neigung von Rutschen, Gleit- 
bahnen usw. sowie des Schiittwinkels von Haufwerken dienen. 
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INTRODUCTION 

Tur angle of friction of dry granular materials 
sliding down chutes, slides, etc., has often been 
measured and reported in the literature, but 
careful perusal of papers reveals no serious 
attempt to investigate the factors that affect 
the angle of friction. The term “ angle of repose ” 
is synonymous with “angle of friction,” but in 
most engineering applications it is used only to 
designate the maximum slope taken up by heaps 
or beds of materials before sliding occurs [1]. 

At the beginning of this century a great deal of 
research was carried out into the effect of earth 
pressure on the angle of repose of soils by WiLson 
[2], Crosstuwarre [3], Breve [4], Darwin [5] 
and Car [6]. The main purpose of this work 
was to investigate the earlier theories of RANKINE 
[7, 8, 9, 10). Rankine does not appear to have 
measured angles of friction himself, as the figures 
he quotes are taken from a previous investigation 
by General Morin who was interested in trench 
embankments, ete., for military purposes. It 
appears, therefore, that the figures for angle of 
friction have been handed down by text books 
for at least a century without a check being 
applied. The table of angles of friction appearing 
in the Handbook of Physics and Chemistry [11] 
of 1953 is compiled from the data of RaNKINE in 
1858, who obtained the figures from other earlier 
investigators. 

Recent work on the angle of friction of sliding 
metal blocks has been carried out by Burwe.. 
and Rapinowicz [12], Parker and Harcu [13] 
and BikerMan [14, 15, 16], who have put forward 
theories 1954 
Ricurer [17] measured the static coeflicient of 


based on Amonton’s Law. In 
friction for chopped grass, corn silage, hay and 
straw 


obtain data for the design of chutes. His method 


sliding on galvanized steel in order to 


of measurement was a variation of the original 
experiment of CouLoms for measuring the vertical 
and tangential forces applied to a frame contain- 
ing the material. Hotsrook and Fraser [18] ob- 
tained some interesting results on angle of friction 
while investigating the sizing of coal by screening. 
They found that the angle of friction for coal 
sliding on bright polished steel was less for round 
particles than for flat or wedge-shaped particles. 


They observed that for larger sizes of coal the 
angle of friction was less than for small sizes, 
and the angle increased rapidly as the size 
decreased and the specific surface increased. 
The object of this investigation is to find the 
effect of particle and surface properties on the 
coetlicient of friction. The properties considered 
are: shape factor, specific gravity and average 
diameter of the granular material , and roughness 
and specific gravity of the surface on which the 
granules slide. It was expected that an empirical 
dimensionless relationship of the form 


Tv tan @ db E . 


would be obtained. 


EXPERIMENTAL PROCEDURE 
RESULTS 


AND 


A wooden chute was constructed of dimensions 
18 in, 
raised at the other by a fine-threaded screw so 


6 in. Gin, pivoted at one end and 


that angles from the horizontal between 0° and 
5° were obtainable. The angles of inclination 
were measured by means of a plumb-bob and 
reading in half degrees but 


allowing estimation to the nearest quarter of a 


graduated = scale 
degree. The floor of the chute was interchange- 
able, so that surfaces of different roughness and 
specific gravity could be used. 

Surfaces were made up by coating a base board 
with glue and sprinkling different materials over 
the glue so that, on drying, a surface of a suitable 
degree of roughness was obtained. 

Different the 
relevant physical properties being determined 
beforehand. 


granular materials were used, 


The average shape factor of the 
particles was determined by permeability methods 
using a modified Lea and Nurse [19] apparatus, 
and by applying Carman’s [20] equation. The 
materials used and their physical properties are 
shown in Table 1. 

The roughness of the constructed surfaces was 
considered to be the average depth of the inter- 
stices between the particles making up the surface. 
This is analogous to the roughness measurement 
in pipes for friction loss determinations where 
the average depth of the pits and irregularities 
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Table 1. 


Physical properties of granular materials 





Shape Average 


Specific gravity at 25 © 
pecific gravity at 25 factor 


s(H,O 1) 


diameter 
D (em) 
00-4026 
O34 
2517 
00-1924 
0-0592 
0-098 1 
0-0661 
00-0508 
O-O3T2 
O-OR2T2 


Material 


11340 
1-380 
144 
1-103 
1-500 
2-647 
2647 
2-647 
2-47 


247 


Lead shot 
Wheat 
Rice 
Rape seed 
Sugar 
Sand A 
Sand B 
Sand © 
Sand ID 
Sand E 





in the surface are taken as a measure of roughness. 
The roughness index (¢) for the surfaces is there- 
fore proportional to the average diameter of the 
particles making up the surface. Ten surfaces 
were prepared of differing roughness and specitic 
gravity. 

Each of the materials of Table 1 were tested 


on the prepared surfaces in an experiment of 


100 runs suitably randomized to prevent bias due 
to wear of the surfaces and experimental error 
from affecting the statistical analysis. 


at which at least 50 per cent of a bed of granular 
materials slide freely down a prepared surface. 
It was found that the bed should be as thin as 
possible so that the maximum influence of the 
roughness of the surface could be exerted. If a 
thick bed was used it was found that the angle 
measured was an angle compounded of the angle 
of repose of the material and the sliding angle. 
The criterion that at least 50 per cent of the bed 
must be discharged at the sliding angle was 
adopted as it was found that some material was 
discharged at angles greater and less than the 
true angle, due to outside influences over which 
there was no control, e.g. accidental vibration 
due to the method of raising the chute. 

The results obtained for angle of friction are 
tabulated in Table 2. Preliminary observation of 
the results in Table 2 
between angle of friction and roughness of surface 


indicated some relation 


for each material tested, also for a given roughness 
of surface the angle of friction tended to increase 
as the particle size decreased. 

It was considered that a relationship of the 


Y 
Ty +6(5) 


form 


p= tand 


The angle of friction was measured as the angle 


Table 2. 


Experimental results for sliding angle for different materials and surfaces 





| 
Roughness Sliding materials 
Sliding 
surface 


(em) 7 
Lead 
Sand D | Sand KB 
41-0 
38-0 


shot Wheat Rice Rape Sugar Sand A Sand Bs Sand C 


Wheat 
Rice 

Rape seed 
Sugar 
Sand A 


39-0 
338-0 
37-5 
38-0 
38-5 
37-0 
37-5 
3o-0 
39-0 


32-5 


41-5 
39-0 
39-0 37-5 
39-0 38-0 
39-0 37-5 
39-0 37-0 
39-5 37° 

38-0 37° 

37-0 37° 

32-0 34-0 


O-3134 | 2240 BB. te 35-0 0 $1. wo 
$1. 


1- 
1: woo 


39-5 


5 
990 | 39-5 
410 39-5 
4-0 39-5 40-0 
42-0 39-0 39-0 
1-5 39-5 39-5 
0-0 38-0 BO-+4 
39-5 38-0 | BR 
31-0 19-5 39-0 37-0 


| 23-0 16-0 30-0 20-0 
| | 


34-5 
33-0 
34:5 


2517 32-5 
0-1924 
00-0592 
0-098 1 
0-0661 
0-0508 
1D O-O372 
» E 0-0272 

Smooth plastic O(approx.) 


| 


3o-5 
3-0 
374) 32-5 
33-0 26-0 
34.0 25-5 


32-0 23-0 


—_ 





Excessive humidity may cause sticking of the materials and consequent alteration of angle of friction. The 
materials were dried and the experiments were carried out with average atmospheric conditions of 70-75 F and 
Due to the influence of the relative himidity, the effects of electrost .tic 


Note 


relative humidity of 50-60 per cent, 
attraction were considered to be negligible. 
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would explain the variations where 2 and y are 
exponents and so the data was analysed and the 
values of the exponents determined. Initial tests 
showed that the intluence of the specitic gravity 
The 


regression analysis was continued and the final 


of the surface was negligible. statistical 


equation derived was of the form 


2110 “ 
p = tan @ = —— - 0-436 (=) 
f.* D 


O-OITL s, | OIS34 (3) 


The equation was found to have a correlation 
coetlicient of 0-92 with a standard deviation of 
0-OS836, 


The value of «, D was assumed to have a maximum 
value of unity as any increase above unity means 
that the particles of the flowing material are 


the When 


the flowing material will fill the interstices and 


smaller than interstices, this occurs 
alter the surface so that the material flowing is 


sliding over a surface of the same characteristics, 





——Proposed equation 

















05 


' Mf $ 


OS 
Pies. 1, 2 and 3. A plot of tan @ against | . ) for wheat rice, rape seed and lead shot. 
D 
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and the angle of friction is equal to the angle of 
repose of the material, 


DiscuSSION AND CONCLUSIONS 


Analysis of the results indicates that the correla- 


tion is highly significant and that the angle of 


friction of the materials considered depends on 
the shape factor, specific gravity and average 
diameter of the flowing particles, and on the 
roughness of the surface upon which the particles 
are sliding. 

The value of « may be considered as a relative 
roughness term and can have values between 
Oand 1. Values of zero are to be expected if the 
surfaces are smooth, as in polished steel, ete., and 
values of unity are to be expected when the 
material flows over a surface composed of materials 
identical to itself, i.e. that the angle of friction is 
equal to the angle of repose, 

The results for wheat, rice, rape seed and lead 


. 
) 


shot are plotted against (« D)"° and are shown in 
Fig. 1, 2, 


fall close to the curve as given by equation (3), 


and 3. The values for wheat and rice 


but the values for rape seed and lead shot are 


lower than predicted due to the probability of 


rolling of the particles rather than sliding. 
The 


angle in accordance with variations as stated by 


equation predicts variations of sliding 


other investigators in the following manner. 


(1) Materials of the same average diameter 


and specific gravity sliding down 


known roughness will have a higher angle of 


friction as the particle shape factor decreases 
from unity, Le. irregular shaped particles slide 
less freely than spherical particles of the same 
average diameter. 

(2) Materials of the same shape factor and 
average diameter sliding down a surface of known 
roughness have angles of friction that decrease 


as the specific gravity increases, ie. heavier 


particles slide easier than lighter particles of the 


same average diameter. 


a surface of 


angle of friction of granular materials 


(3) Materials of same shape factor, average 
diameter and specific gravity sliding down a 
surface of known roughness have angles of 
friction that increase as the roughness of the 
surface increases, i.c. the rougher the surface, the 
vreater the angle for sliding until the relative 
roughness equals unity and the material then 
slides over itself. 

(4) Angles of friction of sliding materials are 
independent of the composition of the surface, 
except for its physical state. It is not possible to 
compare the results obtained in this investigation 
with the results of other investigators due to the 
the results were 


dearth of information on how 


obtained. Rankine (7) publishes angles of 

friction for sand without specifying either its 

properties or method of measurement. 
The this 


indicates the influence of the variables only, and 


equation derived in experiment 
it is not expected that equation (3) will be 
applicable to prediction of angles of friction for 
cases other than materials with the property 
ranges quoted or by other methods of measure- 
ment of the angle of friction. It must also be 
that 


using dry 


noted the experiments were carried out 
and that the 


water can appreciably affect the angle of friction. 


materials presence of 
Alternative methods of measuring the angle of 
friction of materials are being considered and 
tested and it is hoped that, with the new data 
obtained, equation (2) can be rechecked. 


NOTATION 
average particle diameter of granules 
roughness index of surface 
average shape factor of granules 
tan @ coeflicient of friction 


where 0 angle of friction 
specilic gravity of granules 
[relative to H,O 1} 


specitic gravity of surface [relative to H,O 1] 
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Abstract—The effect of feed rate, frequency of vibration, angle of inclination and screen aperture 
on the effectiveness of a single deck vibrating screen was investigated by a statistically planned 
experiment of 256 runs using clean sand. 


It was found that in order to get a high effectiveness it was necessary to select proper combin- 
ations of feed rate and screen aperture as main variables and to control angle of inclination and 
frequency of vibration as the main interaction, but it has been shown that certain combinations 
of these factors are more important than others, It has been found that all things being constant 
the effectiveness of separation increases with : 

(a) an increase in screen aperture 

(b) an increase in angle of inclination up to maximum of about 15° and then progressive 

decrease in effectiveness with higher angles 

(c) a decrease in the feed 

(d) an increase in frequency of vibration to a certain maximum and then a decrease in 

effectiveness with higher frequencies. 
The influence of these factors upon the attrition of the sand particles showed that for sand 
the amount of * fines produced is negligible, indicating no change in the size of the product. 


Résumé—Les auteurs ont étudié les différents effets de la vitesse d'alimentation, la fréquence 
de vibration, langle dinclinaison et le nombre de mailles du tamis, sur leflicacité d'un simple 
tamis vibrant, par une expérience de 256 essais avec du sable propre. Hs ont montré que, pour 
avoir un bon rendement il faut ajuster la vitesse d’alimentation au nombre de mailles du tamis 
comme principales variables et Pangle dinclinaison a la fréquence de vibration, dont interaction 
est importante. Certaines combinaisons de ces facteurs sont plus importantes que d'autres. 
Les auteurs ont montré que tous les paramétres étant constants leflicacité de la séparation croit 
avec : 
(a) une augmentation de Pouverture du tamis. 
(b) une augmentation de Tangle d'inclinaison jusqu’é un maximum denviron 15°, aprés 
quoi le rendement décroit pour des angles plus grands. 
(c) une alimentation plus faible. 
(d) une augmentation de la fréquence de vibration jusqu’é un maximum, car ensuite la 
séparation décroit avec des fréquences plus élevées. 
Linfluence de ces facteurs sur lusure des particules de sable a montré que la quantité de 
“ fines’ produites est négligeable ce qui prouve qu'il n'y a pas de variation de dimension du 
produit. 


Zusammenfassung -Der Einfluss der Aufgabemenge, Schwingungsfrequenz, Neigungswinkel 
und Maschenweite auf die Wirkung cines einfachen Schwingungssiebes wurde durch einen 
statistisch geplanten Versuch von 256 Messungen an reinem Sand untersucht. 


Um einen hohen Wirkungsgrad zu erzielen, musste man geeignete Kombinationen von 
Aufgabemenge und Maschenweite als Hauptvariable auswihlen und Neigungswinkel und 
Schwingungsfrequenz als wichtigste Kinflussgréssen regeIn, aber es wurde gefunden, dass best immte 
Kombinationen dieser Kaktoren wichtiger waren als andere. Leiss man alles iibrige konstant, 
so stieg die Trennwirkung an mit : 





nh. T. 


steigender Maschenweite. 


steigendem Neigungswinkel bis zu cinem Maximum von 15 , 


progressiv abzufallen 


fallender \ufgabemence 


Fownern amd S.C. Lam 


um bei hoheren Winkeln 


(dl) steigender Schwingungsfrequenz bis zu cinem Maximum, am bet hoheren Frequenzen 


wieder abzufallen 


Der Kintluss dieser Grossen auf den 


Abrieh der Sandteilchen zeigte, dass fiir Sand die 


Eerzeugung feinster Teilchen vernachtiissigbar blieb, dass also die Korngrosse des Produkts nicht 


verindert wurde. 


INTRODUCTION 


ScREENING ts one of the most commonly used 
methods for separation of finely divided solids in 
industry. Among the various types of screens 
used, the vibrating screen is the most important 
since 1920) after 


the 


materials to withstand the vibration stresses set 


having made great progress 


successful manufacture of constructional 
up in its operation (1). 

The effectiveness of separation of solid particles 
by a vibrating screen can be affected by many 
different variables, and it is due to the large 
number of these variables that progress so far 
has been restricted in detining optimum operating 
conditions for these screens. The variables that 
can affect the effectiveness of separation are as 
follows, as proposed by many different workers 
[}, 2, 3, 4, 5, 6, 7, 8 and 11). 


(1) Variables duc to the material being screened 
(a) Bulk specific uray ity of feed. 
(b) Particle shape. 
(c) Percentage of near size (0-7—-1.5 times 
the sereen aperture) material in the 
feed. 

(d) Moisture content of the feed. 

(ec) Static charge generation. 

(f) * Stickiness * 

(zg) Abrasion resistance to attrition. 


of the material. 


Mechanical variables due to th 
screen used 

(a) Length and width of screen. 
(b) Amplitude of vibration. 

(c) Frequency of vibration. 

(d) Angle of inclination of screen. 
(e) Direction of vibration. 

(f) Capacity of the screen, ic. feed rate. 


type of 


(3) Variables due to the sercen cloth 


(a) Mesh aperture, 

(b) Free area available in cloth. 

(c) Type of sereen, ie. square, circular, 
rectangular aperture, ete. 

(d) Feed rate to prevent “ blinding.” 

(ce) Material of construction, ie. copper, 

steel ete. to withstand distortion of 

mesh with time. 

Uniformity of mesh size over area of 


screen, 


It is due to this large number of variables that 
for 30 years workers have been prevented from 
evolving a general relationship that would allow the 


prechetion of effectiveness of a sereen. However, 


some general rules have been proposed by several 


workers and these may be summarized as follows : 


(a) A wet, sticky material should be dried 
as much as possible before dry screen- 
ing, although under certain circum- 
stances excess water can be used in 
wet screening. 

The 


vibration depends on the median size 


frequency and amplitude of 
of the particles in the feed and there 
are optimum values that give the 
maximum effectiveness of screening 
[5] [9], but these factors work in 
opposition, 

The angle of inclination of a sereen 
should be such that “ blinding” is 
prevented with too low an angle, but 
time of residence shall be as long as 
possible before discharge [5]. 

The of vibration of the 
screen should be counter-flow to the 


direction 
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material flowing down the screen in 
order to get maximum agitation and 
residence time (5). 

The sereen length has an optimum 
value of about 2-2) times its width. 
Square-aperture screen cloths have a 
greater free area than sereen cloths 
with circular apertures (1). 

The capacity of a sereen, defined 
as lb of feed, ft® of screen area hr, mm 
screen aperture, should be kept con- 
stant for any particular screen to 
prevent “blinding” and loss of 

effectiveness. 

It is the intention of this paper to present 
results of effectiveness of screening of dry sand 
on altering the four variables of feed rate, angle 
of inclination, frequency of vibration and screen 
aperture and to estimate the amount of attrition 
that may be produced as the sand passes through 
the screener. 


EXPERIMENTAL APPARTUS 


The vibrating screen used in the experiments 
was a single deck Denver-Dillon screen vibrated 
by an off-balance flywheel. The speed of vibration 
could be varied by the use of various diameter 
pulley wheels driven from | h.p. electric motor, 
The off-balance flywheel gave a tixed amplitude 
The angle of inclination of the 
to 20. 
over-all size was fixed at 10 in. width and 24 in, 
length giving an effective 
1-76 ft?. The feed rate to the screen was controlled 
by varying the orifice diameter of a bottom 


at all speeds. 
screen could be varied from 0 The sereen 


screening area ol 


discharge conical hopper which had been pre- 
viously calibrated. Sereens of various apertures 
were selected made of copper woven mesh giving 
square openings. 

Clean, dry river sand was used in the tests of a 
size distribution covering the range of screen 
apertures selected so that measurable amounts 


of oversize and product could be obtained. 


EXPERIMENTAL PROCEDURE 

Clean river sand in 5 ton lots was obtained and 
air dried, A screen analysis on this sand showed 
that reasonable proportions of oversize and under- 
size could be obtained if screens of aperture 30, 
40, 50 and 60 mesh were used in the experiments. 
The speed of vibration was fixed at 952, 1130, 1326 
and 1498 rev ‘min by the use of various diameter 
pulleys on the driving motor. The machine had a 
The 
angle of inclination of the screens were fixed at 
6. 11, 15° and 19 
could be 


fixed amplitude of 7 64in, at all speeds. 
while the feed rate to the 
10-0 


5-5, 7-25, 


sereens controlled at 
and 15-45 lb min. 

A statistically 
i + 1 t 


the four levels of each of the four factors were 


planned experiment — of 


256 runs was developed in which 


split arbitrarily into pseudo-factors giving 2° full 
factorial design. This design was split into 4 
blocks of the block 
conpounded with the third order interaction. 


64 runs with differences 
This ensured that any differences in sand quality 


the blocks did 


a bias to the The runs in each 


being used for individual not 


cause results. 
block 
The 


given in Table 1. 


were suitably randomized.” 


values of the factors and levels used are 


Table 1. Table of levels and factors used in the statistical design 





Feed rate 
B 
(Ib) min) 





952 
1130 
1326 
1480 


Speed of 


vibration 


(rev min) 


Factors 


Angle of 


inclination Screen aperture D 


Mesh No. | Aperture (ys) 


60 
30 | 
wo | 


276 
318 
447 


30 506 








K. T. 


Table 3. 
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Analysis of variance for effectiveness 





Nature of 
effeet 


Source Sum of squares 


tre tore a) 


Main factors 
LO4AS-SSU 


sze7-4aze 


22.2 h)0R20 


T25-441 

rior isl 
Interaction ZIGS- 821 
between pairs of 


factors 11887-3501 


16-451 


ihe 23 BSY 
Interaction iBDbp 
between triplets 
of factors ROD 


Cpa 


1206-080 


Residual inCD 3385-173 


Total 57316-348 


DF. 


Variance I ariance 


estimate ratio 


1S O16 o-410* 


Soo 


OS 145 OS25 


7415-343 | 30-7634 


SO-SL iB ABC 0-943 


5117s i ABCD 0-601 


2hos0 1D ARCD 3-766F 


1SVO-S A BC BCD 14053F 


Is BD BCD O200 


CD BCD OTe 


1BC ARCD 2-050T 
BD ABCD 
BCD ABCD 


CDA ABCD 





* significance level of 1 in 20 


Standard error of estimate of effectiveness 


The screener was set to the particular conditions 
required and sand allowed to flow onto the screen 
When 


this occured two samples of both oversize and 


until steady conditions were obtained. 


undersize were collected in 10sec and weighed. 


The samples were reduced in size and then test 


screened on Tyler standard screens in a vibrating 
shaker. The size distribution was then obtained 


for each sample. Since the aperture of the 


+ significance level of 1 in 1000 


(4+ per cent. 


screens in the machine differed slightly from the 
corresponding test sieve the actual mass in the 
undersize was interpolated from the assumption 
that the cumulative mass fraction retained on 
the test screens of apertures above and below 
the actual screen was a linear relationship. This 
assumption was tested and found to be a good 
approximation and obviated the necessity of 


graphically representing the cumulative mass 
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fraction retained against for every 
sample taken. 

The effectiveness [10] of the 
calculated from the relationship defined as the 
product of recovery and rejection of material 


of the desired mesh size. 


aperture 


screen Was 


U 


Recovery 1,0 


(1 
(x, U 


ty) U 


Rejection - ay 0) 


Thus effectiveness °, 


Recovery Rejection 100 (1) 


The mean mesh size of the material leaving the 
screen was calculated as the sum of the mean mesh 
sizes of the undersize and oversize with respect 
to their individual mass flow rates as in equation 
(2). 

M,0O U 
Oo 


M,, 


u (2) 


Table 3. 


The value of M was calculated from the sereen 
analyses of the undersize and oversize by the 
relationship 
_ 2 (mD) 

= (m) 


M-:- (3) 


RESULTS AND CONCLUSIONS 
The complete results for the 256 runs for 
replicate runs are available on request to the 


authors. A statistical analysis of these results 


” 


vielded the analysis of variance given in Table 2. 
The analysis of variance for the total mean 
mesh size expressed in « is given in Table 3. 
Examination of the results for the effect of the 
various factors upon the per cent effectiveness of 
the screens showed that only two main factors 
of screen aperture and feed rate were important, 
with aperture considerably more significant than 
feed rate. Two signiticant first order interactions 
of feed rate with aperture and frequency of 
vibration and angle of inclination were obtained, 


Analysis of variance for mean mesh size 





Nature 


of effect 


Source 


Main 
factors 


AB 
AC 
AD 
BC 
BD 
CD 


Interaction 
between pairs 
of factors 


AABC 
BD 
BCD 
CDA 


Interaction 
between triplets 
of factors 

Residual ABCD 


Total 


Sum of 


Muares 


10.358 


} artance 


estimate 


730 
376 
zt 
1459 


$556 
2511 
1679 
this 
3428 


2251 


te t 


a) =) oh 


TSS0 
780 
3429 


tt 


$569 


23 








No factor or interaction proved to be significant. 


Standard error of estimate of mean mesh size — Jl 
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Fic. 1. 
effectiveness of a screen with increasing level of the main 


Graph showing the trend of the average per cent 


variable. 


while the only significant second order interactions 
were feed rate, frequency of vibration with angle 
of inclination and frequency of vibration, angle 
of inclination and aperture. Since the design of 
the experiment was such that the third order 
ABCD with the 


residual, the above conclusions must be regarded 


interaction was confounded 
as relative as it may well be shown that a combina- 
tion of all four factors could affect the effectiveness 
of the screen. 

A clearer interpretation of the results can be 
seen from Fig. 1, 2 and 3. Fig. 1 shows the effect 
of the four factors with increasing level on the 


e130 
s ° ¢ 


Fig. 2. 
and frequency of vibration on the per cent average 
effectiveness. 


Graph showing the effect of angle of inclination 


average per cent effectiveness if operating con- 
ditions are regarded as constant. Fig. 2 shows 
the effect on the effectiveness with an alteration 
in both feed rate and aperture while Fig. 3 shows 
the effect of angle of inclination and frequency of 
vibration on the effectiveness. 


Fic. 3. 
of screen on the per cent average effectiveness. 


Graph showing the effect of feed rate and aperture 


As a check to the findings of the main experi- 
ment and to estimate the standard error, a set 
of duplicate runs was performed. The 9 runs were 


performed on 3 levels of frequency of vibration 


and angle of inclination, with the feed rate and 
fixed at 7-25lb min and 30 
The analysis of variance for both 


aperture mesh 
respectively. 
effectiveness and mean mesh are shown in Tables 
4 and 5. 

The results of Tables 4 the 
conclusions obtained from Tables 2 and 3. The 
estimate of the standard error from Table 2 is 
0-4 per cent for the per cent effectiveness which is 
higher than that of Table 4 of 0-24 per cent, 
but this may be due to the fact that Table 2 is 
based on 256 runs while Table 4 is based on 18 
and block differences are included in the residual 
for Table 2. The estimate of the standard error 
in mean size for Table 3 and Table 5 are 1-1 « and 
1-8 respectively. This is indicative of the fact 
that the amount of attrition is negligible and the 
change in mean mesh size of the material passing 
This conclusion 


and 5 confirm 


the screen is not significant. 
only applies to the sand used in the experiment 
and it is obvious that this will not be true for all 
materials, 
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Table 4. Analysis of variance for effectiveness for 18 runs 





Sum of 
squares 


Nature of 
effect 


Source 


237-23 
Main factors 





an 


Interaction | BC 
between B and ¢ 


Replications Residual 9.01 


Total 383-45 


Variance Variance 


estimate ratio 


118-62 


21-29 





*Significance level 1 in 1,000. 


Standard error of estimates of effectiveness = 1-00/ 4/18 


Table 5. Analysis of variance 


== (24 per cent. 


for mean mesh size for 18 runs 





Sum of 


sy mares 


Nature of 
effect 


Source 


Main factors 


Interaction 
between B and C 


Residual 


Replications 


Total i710 


Variance 
ratio 


Variance 
DF. 


estimate 


Residual 
c 
Residual 





No factor or interaction proved to be significant, 
Standard error of estimate of mean mesh size 


In the runs performed the percentage effective- 
ness of the screen using sand varied from a high 
value of 94.5 per cent for a feed rate of 7-25 lb min, 
952 rev min, 15° angle of inclination and an 
aperture of 30 mesh (run No, 231) to a low value 
of 22-5 per cent with identical conditions except 
that the angle of inclination was 6 degrees 
(run No, 67). This is probably due to * blinding ” 
of the sereen as at low angles the oversize material 
cannot escape quickly enough from the screen 
surface. Examination of Fig. 1 shows that, for 
all conditions being kept constant, the 
cent effectiveness of a screen increases as: 


per 


614/18 = 18 pw. 


(a) the feed rate decreases 

(b) the angle of inclination increases to a certain 
value of the angle when the effectiveness 
decreases with further increase in angle 

(c) the frequency of vibration increases and 
then decreases as the frequency increases 

(d) the screen aperture increases. 


The main conclusions developed from this 
paper were used as a basis for a set of new experi- 
ments using different materials and varying the 
amplitude of vibration. It is hoped to present 
these results in a later paper. 
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NOTATION 


aperture of screen which just prevents passage of material 


feed rate of material 

mass of material retained on screen of aperture D 
mean mesh size of material 

mean mesh size of oversize material leaving sereen 
mean mesh size of undersize material leaving screen 
oversize rate of material 

undersize rate of material 

mass fraction of desired mesh material in oversize 


mass fraction of desired mesh material in undersize 
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Research on ion exchangers by means of the percolation technique 
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Abstract—In ion exchange important process characteristics are the equilibrium curve, relating 
the concentrations in the two phases, when they are in equilibrium, and the kinetic data, governing 
the mass transfer, when both phases are not in equilibrium. It is shown, that these data can be 
determined by simple percolation experiments, 

The first part of this paper describes a method to determine the equilibrium relation graphically, 
that is generally valid if the equilibrium curve shows no inflexion point. For exchange of ions 
of equal valency a numerical method is deduced, which also gives the bed capacity. These methods 
are applied to a number of systems. It is shown that the anion species exerts a strong influence 
on the exchange if hydrogen is one of the exchanging species. This phenomenon is explained semi- 
quantitatively. 

The second part describes a method to determine the mass transfer characteristics graphically 
from a sharp break-through curve and the equilibrium curve. The break-through curve is of 
the stationary type, in which case the relative concentrations in the two phases on one cross- 
section are equal, 


This method is applied to a number of systems with varying flow rate and particle diameter. 


Résumé Les caractéristiques importantes des processus d’échanges ioniques sont la courbe 
deéquilibre reliant les concentrations dans les deux phases ct les données cinétiques qui gouvernent 
le transfert de masse quand les deux phases ne sont pas en équilibre. Ces donneées sont déterminées 
dans une simple expérience de percolation, Dans une premiére partie, Pauteur donne une méthode 
graphique de détermination de courbe d'équilibre qui est généralement valable lorsque la courbe 
ne présente pas de point d'inflexion. Pour les échanges d’ions de méme valence Pauteur déduit 
une méthode numérique qui donne la capacité du lit. Ces méthodes ont été appliquées a un 
certain nombre de systémes. Les anions exercent une grande influence sur léchange surtout 
en présence d'ions hydrogénes. Ce phénomeéne a été expliqué semi quantitativement. 

Dans une 2e partie lauteur décrit une méthode de détermination graphique du transfert de 
masse & partir de la courbe d’équilibre et d'une courbe nettement brisée. La courbe brisée est 
du type stationnaire, et dans ce cas les concentrations relatives dans les deux phases pour une 
section sont égales. 

Cette méthode est appliquée a un certain nombre de systémes en utilisant des écoulements 


et des diamétres particules variables, 


Zusammenfassung —Wichtige Kinflussgrossen fiir den lonenaustausch sind die Gleichgewichts- 
kurve, die die Konzentrationen in den beiden Phasen im Gleichgewicht enthalt, und die 
kinetischen Daten, die den Stoffiibergang im Ungleichgewicht wiedergeben. Es wird gezeigt, 
dass man diese Werte aus cinfachen Perkolationsversuchen erhalten kann. 

Der erste Teil dieser Arbeit beschreibt cine Methode zur graphischen Bestimmung der 
Gleichgewichtsbeziehung, die allegemein giiltig ist, wenn die Gleichgewichtskurve keinen 
Biegepunkt aufweist. Fiir den Austausch von lonen gleicher Valenz ist cine numerische 
Methode abgeleitet, die auch die Bettkapazitét ergibt. Diese Methoden werden auf eine grosse 
Zahl von Systemen angewandt. Es wird gezeigt, dass die Anionenart einen grossen Einfluss auf 
den Austausch ausiibt, wenn Wasserstoff eine der austauschenden Arten ist. Diese Erscheinung 
wird halb quantitativ erklirt. 


*Present address: Dorther Weg 16, Epse, Near Gorssel, Netherlands, 
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Der zweite Teil beschreiht cine Methode zur graphischen Bestimmung der Stoffiiber- 


gangsygrossen mit Hilfe einer scharfen Durchbruchskurve und der Gleichgewichtskurve. 


Diese 


Durchbruchskurve ist stationdir, in welchem Falle die relativen Konzentrationen in den beiden 


Vhasen auf cinem Querschnitt gleich sind. 


Diese Methode wird auf cine grosse Zahl von Systemen mit veriinderlichem Durchfluss und 


Teithchendurchmesser angewendet 


l. INTRODUCTION 


In ion exchange, as in all mass transfer processes, 


both the equilibrium and the kinetic character- 
istics of the process are of interest. As a conse- 
quence, it ts important to be able to determine 
these characteristics in a simple way. 

The 


processes is mostly determined by means of batch 


equilibrium relation of ion exchange 


experiments. A modification of this method 
exists in which a column technique ts used. Here 
a column of ion exchange particles is equilibrated 
solution inp which the 


and B 


By regenerating the column contents one can 


(by percolation) with a 
concentration ratio of tons A is known. 


determine the concentration ratio of these tons 
in the resin phase, Then one point of the equilib- 
rium curve of ion pair A* — Be is known. 

The kinetic characteristics of an ion exchange 
process are usually determined by measuring the 
slope of the break-through curve resulting from a 
percolation expermment, in which the displacing 
ions show greater allinity towards the resin phase 
than the ions, originally bound by the ton ex 
changer. Although this method is used almost 
exclusively, there exists some confusion im the 
interpretation of the measured break-through 
curves. 

In this paper it will be shown, that for mono- 
functional the 


relation can be determined from one break-through 


exchangers binary equilibrium 


curve. Further the determination of the mass 
transfer coetlicients from the equilibrium relation 
and a break-through curve will be treated. 


2. Turory 
The following assumptions are made concerning 
- 
the packed bed and the flow through it. 
(a) The ion exchange particles are supposed to be 
uniformly packed ; 
(b) Swelling and shrinking of the particles is 
neglected at present ; 


(c) Longitudinal mixing is supposed to be neg- 
ligible. 


With these assumptions a material balance 
around an infinitesimal slab of the bed gives the 


equation ; 


Introducing x = C,C, and y = q qd, the relative 
the liquid and resin phase 
R F, the superticial flow 


velocity, equation (1) becomes : 


concentrations mn 


respectively, and v, 


ir F Jo OY E 


€ = 0 
” Ug dt 


Instead of the independent variable t one can 
introduce ¢ — (€ = vu); the term ez v, is the time 
to fill the the 
particles with the solution. So 
[tf — (ez v)| Col is the total number of equivalents 


that has flown out of the column since t 0. 


necessary voids between resin 


Incoming 


After the following substitutions : 


T ( . =} C,R, and 
“oe 


s=q Fz total bed capacity, 


there results: 
\ ‘ 
(5) (<") -0 (3) 
os), dt]. 


When for ¢ <— 0 the bed is completely charged 
with B* ions and washed with pure water, and if 
the incoming solution contains A* ions only, in 
concentration Cy the boundary conditions for 


equation (2) are: 
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EquiLiprium 
EXPERIMENTS 


3. DETERMINATION OF THE 
RELATION FROM PERCOLATION 
GENERAL SOLUTION 


When in an ion exchange reaction only two ion 
species, A* and B* are present, then 1 and y are a 
direct measure of the composition of the two 
phases. 

It is assumed that the equilibrium relation is 
given by: 


y = f(a), (5) 


and that £" (x) ( da*) is nowhere equal to 
zero. 


The latter is generally the case with mono- 


df (x) 


functional exchangers. 

It is well known from literature, [8, 18], that 
under suitable conditions the experimental break- 
through curve is nearly indiscernible from an 
“equilibrium ” break-through curve, that would 
result from a percolation experiment, in which the 
two phases on a cross-section are in equilibrium, 

It can be shown [20], that this situation applies, 
if (for f" (27) > 0): 


=f" (2) K,S 


vy f (a) (6) 


The condition f (27) > 0 means that the dis 


placed ions have a stronger allinity towards the 
It can be 


resin phase than the displacing tons, 
shown that the actual and the “ equilibrium ” 
break-through curves differ most for very low and 
very high values of x. 

This criterion (6) and another one (39), to be 
given later on, are easily deduced from the material 
balance (2). 

Fora complete descript ion of this deduction the 
reader is referred to {20}. 

When the criterion (6) is fulfilled it is allowed to 
substitute the equilibrium relation (5) in’ the 
material balance (3). This gives: 


(~  f' (a) (~") 0 (7) 
os], dt] s 

Now a new variable is introduced, 1 — r/s; 1 is 
the dimensionless ratio of the number of equiva- 
lents, that has passed through the column and the 
total bed capacity. 


Introduction of 1 in (7) results in: 


ot’ (2) l}=0 (8) 


This formula has already been given by Siuutn 
[17]. 

In (8) dx/dl is the slope of the break-through 
curve inanw—Il graph. Since dx/dl ¢ 0, we have 
from (8) 

f" (a) (9) 
f (x) Ida (10) 


and 
The meaning of this equation is clear from Fig.1, 


where the hatched area corresponds to y, f (7). 


i. 1. Determination of f (2) from an * equilibrium ” 


break-through curve. 


Thus equation (10) provides a simple way to 
determine the equilibrium relation, either by 
numerical or the 
* equilibrium ” break-through curve according to 
Fig. 1. 

In using this method it is necessary to know the 


mechanical integration of 


total bed capacity, in order to be able to determine 


the value of 


[t (ez Uy) Cy v| 
total bed capacity 


In principle it is possible to determine the bed 
capacity from the same break-through curve by 
integrating it froma — Otox = 1. A much more 
accurate determination is based upon the sharp 
break-through curve resulting from a percolation 
experiment in which the displacing ions show 
greater aflinity to the exchanger than the dis- 
placed ions (f" (a) < 0). 
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This 
f"’ (x) 


equal valency exchange. 


method is generally valid (provided 


> 0), also for reactions where ions of in- 


Numerical solution 
If two ion species of equal valency exchange on 


a monofunctional exchanger, it is possible to 


describe the reaction : 


Ai + Bt At + Bt (11) 


with an equilibrium constant A, defined by : 


[Ap] [Bi] 
(A, ] [Ba] 


brackets 


(12) 


the indicate the 


activities. If. as an approximation, concentrations 


Here square ion 


are used instead of activities, (12) reduces to: 


y(1 r) 
r(l 4) 


(13) 


where » and y are the relative concentrations of 


the displacing tons. 
In this case the general equilibrium relation (5) 


takes the following shape : 
f(r) 


A, 


f' (x) : 
(Ay l ry 


and 
Substitution of (15) in (9) gives: 
(16) 


This equation was also arrived at by Hresrer and 

VERMEULEN {10}. By means of equation (16) it Is 

possible to determine the equilibrium constant 

K,, as well as the bed capacity numerically from 

one break-through curve, for which f (r) > 0. 
Equation (16) is quadriatic in y A, : 


d(ry Ky. d) 


Aywyl \ A, - (1 ryvyl Oo, (17) 


ta (1 r) 
27 \ l 


with = 4/R, Eys (18) 


If the discriminant in (18) is put equal to zero, 


it is clear that d — 0 leads to the condition: 


d¢ 
4 ‘A, r,t 


From ¢ 


o¢ 
“4 " (;*), ai” 


db dM 
of. — d A ' dl 0 
) is K.) \ ; (3) mv hy 


For constant /, (20) can be written as: 


Md dd or 
(se) <4 MR Fry, I 


db 
( \ x) 


(19) 


0 it follows that: 


(20) 


al 


Combination of (21) and results in 


(; UR), - 


wr 
| mane 0. 


(22) 


This equation locates a minimum of 4 aS a 
function of A, with constant 1. This means that 
for each value of 1 (> 1) there is one value of A, 
for which « is smaller than for any other value 
K,; for another value of / there is also another 
value of A, for which x is a minimum, 

Apparently there exists an envelope that is 
tangent to all equilibrium break-through curves 
l graph. 

The envelope is given by the equation : 


in ana 


1 
(23) 
tr (1 — 2) 
An experimental break-through curve gives a as 
a function of the volume of solution that has 


passed through the column, V; in order to cal- 


en 
culate A, from (18) one should know 2 as a function 


of 1. 


found. 


So a relation between V.. and 1 has to be 


er 


The following symbols are introduced here : 


the value of 1, where the equilibrium curve 
tangent to the there 


vj. 


where / 


is a envelope ; 


l= [4r(1 


value of V 


er 
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. 


os 


- on the “ equilibrium ” break-through curve. 


So V,. is the factor that links VU, and / by the 
equation : 


(24) 


Now A, and A, are real, so from (18) : 


l tla , 0 or 


for! 4 1, 
tr (1 


1 
tr (1 vr) 


for] 


Combination of (24) and (25) gives: 
ye 4a(l —2) < Ve forl ¥ l,, 
(26) 
Vi 4a(1 —2) = VW forl ¥ 1, 
From (26) it is apparent that the function 
Ve 4ar(l 


measurements, shows a maximum, and the value 


xv), which can be calculated from the 


of that maximum is Vj... The bed capacity is V.- C5. 

It is then possible to calculate a value of K, for 
every measured point of the break-through curve. 

Relation (18) gives two values for 4/A, (and 
for K,). The highest value should be taken for 
1 >, the smallest one for |< l,,. 

Another procedure to determine the equilibrium 
relation from a break-through curve, for which 
f"’ (2) > 0, was described and used by GLUECKAUF 
[8] and by Duncan and Lister [5]. 


4. Resutts EQuiLipriumM MEASUREMENTS 


In this section some equilibrium data are pre- 


OF 


17 


) 


sented which were obtained with methods des- 
cribed earlier in this paper. 
The 


equilibrium 


the 
sufliciently 


break-through curves from which 


data were calculated 
approached condition (6) of * equilibrium” be- 
tween the two phases. The reader is referred to 
the end of this paper for some remarks on the 


experimental technique used. 


- Nat (Cl-) 


To show the influence of the fluid flow rate on 


The system He Dowex 50 8. 
the shape of the break-through curve, Fig. 2 gives 
the curves for three different rates. It is evident 
that the deviations from the equilibrium break- 
through curve are the most serious in the beginning 
and in a smaller degree in the end of the curve. 
In Fig. 3 the equilibrium curve of this system is 
given in a yw diagram. A, was calculated with 


4 
A 4-00675 cm/sec 
=0133 on/sec 
*%=0°210 cm/sec 
+ %=0-0868 crn/sec 
¥%=0-051 cm/sec | 


“CS O4 —3s — Os 
u 


4 


Equilibrium curve of the system Na‘ 
Dowex 50 8. 


Fig. 3. Hi’ (cl) 


« 
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the method, described in the foregoing paragraph ; 
y was determined from (13). 


The value of K,, calculated from the equilibrium 
1; the equilib- 
> 1; 
therefore the equilibrium diagrams determined 


break-through curve is always - 
rium curves are normally given for A, 
using the theory given above, are turned over an 
angle of 180°. 
gram of Fig. 3 and with all the equilibrium dia- 


This has been done with the dia- 


grams to come. 
From this graph the following conclusions can 


be drawn : 


(a) The shape of the equilibrium curve is in- 
fluenced only slightly by the flow rate. 


observation points ts 


(b) The spread of the 
acceptably small except in the vicinity of 
r — 0-4; the latter phenomenon is caused by 
the fact, that the discriminant of (17) is zero 
there. As a consequence, small experimental 
errors produce large variations in the equilib- 
rium constant. 
The system H* — Na* (SO?) — Dowex 50 8 
The equilibrium curve is again calculated from 
the * 


pared with the curve, determined from batch 


equilibrium ” break-through curve and com- 


experiments, in which a known quantity of the 
exchanger is equilibrated with a known quantity 
of solution. Fig. 4 gives the results for this 
system. 

From this graph the following conclusions can 


be drawn : 

(a) The percolation experiments give equilibrium 
curves which coincide very well ; the tlow rate 
was varied a factor 2 and the experiments 

were performed by three different observers. 

The 

which means that A, is not a true constant 


equilibrium curve is asymmetrical, 


but depends on 2. The question arises, 
whether it then is allowed to use the numerical 
method to calculate A,; in the derivation of 
this method A, was assumed to be constant. 
Comparison with the results of the batch 
experiments shows, that for the (small) var- 
iations in A,, encountered here, the numerical 
method gives results which are equivalent to 


WereveErs 

those of batch experiments. As a conse- 
quence, it is concluded that the numerical 
method may be used for a variable A,, if the 
variations are not excessive. 


The equilibrium curve for the SO,” system 


lies at a greater distance from the diagonal 
y=2 the Cl- curve (Fig. 3). This 
difference is caused by the fact that HCI is a 
stronger acid than H,SO,. This effect will be 
discussed more quantitatively below, 


than 


ae 
sit 


va 





Fic. 4. Equilibrium curve of the system Na* 


Dowex 30 & &, 


Hr (SO?) 


The system H* — Ag* (NO,~) — Dowex 50 x 8 
In Fig. 5 the equilibrium curve of the above 
In the three 


the concentration and the flow rate were slightly 


mentioned system is given. runs 
varied, Further it is apparent that this curve is 
also asymmetrical. 


The influence of the anion on cation exchange 


If the anion present in the solution (Z~) forms 


a weak acid with H* ions, then two reactions 


occur at the same time: 
M; H, == MR, 
H, + Z, ~.HZ 


Hy, 
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Fic. 5. Equilibrium curve of the system Ag 


Dowex 50 » 8. 


H+ (NO,>) 


If HZ is weakly dissociated, reaction (28) will 
withdraw H* ions from the solution and reaction 
That 
means that in this case the charging of the ton 


(27) will be drawn to the right-hand side. 


exchanger with M* ions will proceed further than 
if HZ is completely dissociated, 

This effect will be demonstrated quantitatively 
in the first place for the exchange of Na’ 
H* ions 50 x 8 with either Cl 
SO,” anions, 


on Dowex or 
It is assumed, that there occurs complete dis- 


of HCl of the H* ion 
sulphuric acid. For the dissociation reaction of 


sociation and first of 


the second H* ion, 
HSO, => H* = so, 
the equilibrium constant is given by [15] 


H+} (so? 
i (H+) | Val == 210-2 kg mol m® 


- 29 
. (HSO,) (29) 


After the introduction of concentrations instead 
of activities, (29) becomes : 


[H+] [SO.?) 
(HSO,) 


2:10°-2 kg mol/m® = (30) 


in which y,4 is the activity coefficient of ion A*. 
The activity coetlicients of the univalent ions are 


1 


taken to be equal and ySO¥ is found by extra- 
polating the values given by Lewis and RanpaL 
[14] to higher ionic strengths corresponding to 
0-4 N solutions; in so doing one finds y so,> » 0-2. 
So 


[H+] [SO,* ] 
[Hso,) 


~ 0-1 kg mol ‘m$ (31) 
With (31) it is possible to calculate the true value 
of [H*] instead of the value [H* HSO,-] 
found by titration. Then it is possible to make 
a graph of y against the true value of « 

Na; (Na; + H7) instead of the fictitious value 
Naj (Naz +Hf+HSO,_ ) which is directly based 
on the titration. The result of thiscorrection is given 
in Fig. 6. It is seen that the corrected SO,? curve 
coincides very well with the measured Cl~ curve. 
In the same Fig. the equilibrium curves for 
Na* —H* exchange are given for the 
CH,CICOO™ and CH, COO~, For the correction 
of these curves the following values of the disso- 


anions 


ciation constants were taken [15]: 


0-* kg mol, m$ 
0-* kg mol, m? 


Monochloracetic acid A, 55. 
Acetic acid ; . 





Fic. 6, Equilibrium curve of the system Na* — H* — Dowex 


50 & 8 with different anions. 


The author could not find data concerning the 
activity coeflicients of the anions used here. So 
it was assumed that the influence of the activity 
coetlicients were the same in both cases. 

Then this influence was chosen at such a value 
that the curve for CH,Cl COO coincided approx- 
imately with the Cl- curve and then the same 
value was used for correcting the CH,COO~ curve. 
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It is seen that this last correction fits the CH,COO 
curve reasonably well to that for Cl-. 

Although the latter procedure is not exact, the 
following conclusions can be drawn here : 


(a) In cation exchange the type of the anion can 
exert a strong influence upon the equilibrium 
curve, if H* is one of the exchanging ion 
species. 

It is possible to explain and predict this in- 
fluence by taking the dissociation of the acid 
into consideration. 


These conclusions are a direct contradiction of 


those of Dranorr and Lapipus [4]. 


The exchange of divalent against univalent ions on 
Dowex 50 X 8 
The exchange of reaction of a divalent cation 
against for instance H* ions 


Me}. 2H, =" Me} 2H, (32) 


can, analogous to (12) and (13), be described with 
an equilibrium constant K,, defined by : 


[MeZ"] (Hi? 
(Me? *}[Hy }? 


or 


ry? 
" (3333) 
y)* 


Coy (l 
Go tl 


It is impossible to solve (9) analytically with 
the equilibrium relation (33) as has been done for 
uni-univalent exchange. Therefore Huesrer and 
VERMEULEN [9] try to describe reaction (32) with 
an equilibrium relation of the form of (13). If 
this would be allowed, one could again calculate 
the equilibrium relation from a_ break-through 
curve. 

If, on the other hand, such a procedure is not 
allowed, the equilibrium curve can be determined 
by (for instance mechanical) integration, as pointed 
out in Fig. 1. 

A third method to obtain the equilibrium curve, 
is provided by batch experiments. 

These three methods have been applied to the 
system H* - Zn?* (SO? )- Dowex 50 © 8. 

The results are shown in Fig. 7. 

The drawn line in Fig. 7 is a graphical represen- 
tation of the equation : 


oy “A “+ 2. 


Fic. 7. Equilibrium curve of the system Zn?* — He (SO.2>) 


Dowex JO © 8, 


From Fig. 7 
drawn : 


the following conclusions can be 


(a) It is not allowed, to describe the exchange of 
divalent against univalent ions by an equilib- 
rium constant according to (13). 

The equilibrium curve found by the graphical 
procedure agrees well with the results from 
the batch experiments. 

If this equilibrium is described by relation (33) 
there results a reasonably constant A,. In 

2-1 kg 


this case (Cy — 0-392 kg eq. m*; qo 


eq. m® bed volume) Ky & 7-5. 


Some other di-univalent systems on Dowex 50 * 8 
Of three other systems the equilibrium curves 
were determined graphically, Those systems were : 
H* — Zn®* (Cl-) — Dowex 50% 8; H*~— Ba** (CI-) 
Dowex 50 «8 H+ — Me** (Cl-)~ Dowex 
50 « 8. The Mg system was measured at two 
different concentrations. 
The results are given in Figs. 8, 9 and 10. The 


and 


lines drawn in these diagrams correspond to the 
following values of Ky. 


H* — Zn** (Cl-) - Dowex 50 & 8: 
: B28 at Cy — 0-422 kg eq. 
Dowex 50 « 8: 
19-6 at Cy — 0-592 kg eq. 
Dowex 50 « 8: 
2-00 at Cy — 0-386 kg. eq. 
3-50 at Cy — 0-656 kg. eq. 
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From the Figs. 8, 9 and 10 it is apparent, that the 
value of A, defined by (33) is independent of x for 
the systems studied here but slightly dependent 
on the total concentration. 


7 


ae) 
x 


Fic. 8. Equilibrium curve of the system Zn? HE (Cl) 


Dowex 50 & 8. 


Fic. 9 Equilibrium curve of the system Ba®+— Il’ (Cl) 


~ Dowex 50 8. 


Equilibrium curve of the system Mg?’ ~ H’ 
(Cl-) — Dowex 50 x 8. 


BETWEEN KINETIC 
SHAPE THE 
Break-TurouGcu CurRVES 


RELATIONSHIP 
PROPERTIES 
SHARP 


AND OF 


In paragraph 2 a material balance (2, 3) was 
deduced and the boundary conditions (4) were 
given, 

To completely describe the percolation process 
one further equation is needed. This relation is 
provided by the kinetic equation, giving the rate 
of mass transfer, if there exists no equilibrium 
between the two phases. 

GILLILAND and Bappour [6] and Hrestrer and 
VerMeULEN [9] have used the following equation ; 


kC, ly -y)ax 2 a 


(37) 
| K, ; 


ny 


It is analogous to rate equations of chemical 
reactions, and as such it is not very satisfactory. 
It is to be expected, and also borne out by ex- 
periments, that the rate of exchange is not deter- 
mined by the very fast ionic reactions, but by 
diffusion phenomena in both the liquid and the 
resin phase. On the other hand (37) is the only 
rate equation for which it is possible to solve the 
set of equations (3) and (37) exactly [9]. 

It can be shown that the asymptotic solution 
(e.g. for infinite rate of mass transfer, or infinite 
bed length) for A, 
through curve. 


- 1 gives a symmetric break- 
In view of the fact that such a 
break-through curve is more the exception than 
the rule, one may con-lude, that (37) does not 
describe the mass transfer adequately. 

Another rate relation, based on the diffusional 
character of the mass transfer, is the following: 
or 1 (™ 

Ss 


dM hy, (¢ 


Ci) = ke (Gi —9) 


or 


go —% dy 


S\ot hy, Colv — ai) 


knGo(Y¥i —Y) 

mass transfer per unit area (kg aq m®* sec), 
specific surface of the ion exchange par- 
ticles in the bed (m?, m*). 


If the bed consists of equal spheres S = 6 (1 — e) ‘d. 


k, — mass transfer coeflicient in the liquid phase 
(m, sec). 





mass transfer coeflicient in the resin phase 


Kp 
(mi sec). 

In (38) C and qg (or x and y) are the mean con- 
centrations in the liquid or resin phase. The index 
i indicates the situation at the interface between 
the two phases. 

There are some objections against the use of 
relation (38) too. This type of relation is mostly 
and very successfully used for those heat or mass 
transfer processes, where the temperature or the 
thin 
boundary layer adjacent to the interface and be- 
This situ- 


concentration changes very rapidly in a 


comes constant at a greater distance. 
ation does not occur in an ion exchange process, 

In the resin phase the exact description of the 
unsteady state diffusion phenomena gives a 
solution which is totally intractable [1}. 

In the liquid phase the flow rates used here are 
so low that the flow through the bed ts essentially 
laminar; consequently, the concentration will 
change gradually in a direction perpendicular to 
the direction of flow. 

However, the fact that the use of (38) results in 
consistent values of the coetlicients ky, and k, 
makes it attractive for practical purposes. 

Giuveckaur [7] has proposed other rate re- 
lations ; these empirical equations will not be 
applied here because they lack any theoretical 
back-ground and—being non-linear —are more 
difficult to use, 

Hresrex and Vermevien [9] have tried to 
adapt the exact solution, valid for relation (37), 
to relation (38). There are two objections against 
this procedure ; (a) it is not exact ; (b) it is only 
relatively simple if the mass transfer resistance lies 
mainly in one of the two phases. If both phases 
comparable resistances their procedure 
becomes very cumbersome. 

To be able to solve the set of equations (3) and 
(38) Kuamer [10], [11] and [12] proposed to use 
the stationary sharp front ; that occurs (for A, > 1) 


for: 


show 


qo/ Uo = (8 /M), (39) 
Ve 

so for high values of the residence time v, = [20]. 

This front is characterized by the fact that its 

shape does not change any more, once it has been 


Wevers 


established. This means that the speed, with which 
the front passes through the column ()z/)%),, 
is independent of the concentration x. This speed 


is given by: 


Yo (40) 


M € + qo Co(y/™), 


r 


which is easily deduced from (2). 
The front is stationary if (dv,/dr) 

(dy &), 

From (dy, \”), 


0 or if 
constant. 
C, follows : 


ye Crt, 


Insertion of the boundary values 2 = 0 ; 
r l 
results in: 


(42) 


y r 


So if a stationary sharp front exists, the relative 
concentrations in both phases of a cross-section 
are equal. A similar relation has been deduced by 
Kiamer [10], (11). 

Combination of (42) with (3) provides another 
way of writing this condition for the occurrence of 


this type of front : 


| = 
d x r 


By means of this relation, it is easy to show that 


(43) 


the asymptotic form of the analytical solution, 
given by Hiestex and VermMeuten [9]; 


+4 
1 + exp [k/K,R(I (4*) 


K,) (+ 


s)] 


describes the stationary sharp break-through 


eurve. Also Gruuitanp and Bappour [6] use 
essentially this type of front. 

When the stationary front is used, the relations 
describing the mass transfer, are greatly simplified 


by equation (42). From (38) and (42) follows : 


(45) 


| _ ) 7 k, S(a —.2,) 


where (dv 4). is the slope of the break-through 


ky S (yu, — 9) 


curve, 
It is easily seen that (45) is still further simplified 
if the mass transfer resistance lies mainly in one 


of the two phases. 
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If the liquid phase resistance predominates 
) 
dt). 


f (w *)) and if the resin phase re- 


OK S (a r*) 
Go 


(46) 


[v* from y 
sistance prevails 


7] (47) 


or 
( ) ky S {t(a) — y} 


In Figs. (11) and (12) it is clear that AA’ re- 
presents the driving force in (46) and (47) ; further 
it is apparent that the inflexion point in the break- 
through curve will have concentration 2, (see Figs. 
11 and 12) if the tangent in B is parallel to the 
diagonal. 


[| / 
° et 


‘ 
x 


Fic. 11. Resistance in the resin phase. 


Fic. 12. Resistance in the liquid phase. 


The procedure given above has been described 
earlier by Kuamer [10}. 

If the resistance is in both phases of comparable 
magnitude, it is advantageous to introduce the 


slope m of the line connecting the origin with the 
point A (7;, y;,) which indicates the concentrations 
at the interface. So m = y;/2;. 

Introduction of m in (88) results in ; 


og” Yo (24) ” CoK,(2 —* y 


S \ot m 
Go Ky (ma —y) (48) 
= | e+ 
(K, =(- +—*- 
bien. ky = 
over-all coetlicient 


mass transfer 


based on the liquid phase (m ‘sec). 


Mo Dy 
Tr 
Chk, kr | 


= over-all 


and Ay 


mass transfer coetlicient 


based on the resin phase (m/sec). 
I 


In Fig. (13) it is seen that the driving forces 


x —(y/m) 
Further the ratio of the 


occurring in (48) are indicated by BC 
and by DB = mr — y. 
resistances R, and R, in the two phases 

Ry, =m Kn > 


Rp 
is given by EB/CE 


m cotg o 


(49) 
hey 0 


FD/BF in Fig. (13). 


Fic. 13. Resistance in both phases. 


From the foregoing one may conclude : 


(a) the concentration x, of the inflexion point of 
stationary sharp break-through curve, gives 
an indication what phase causes the biggest 
mass transfer resistance ; 
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(b) since in (49) 4 is nearly independent of # and 
the 
changes when during an experiment & varies 
from 0 to 1, 


m varies with e. ratio of resistances 


If the equilibrium curve can be represented by 
0, to \ kK, 


So. for a large 


relation (13), m varies from A, for 4 
for’ = 1 — yand to 1 fora 1. 

value of A, (both ions differ greatly in their 
affinity towards the resin) it is quite well possible 
that at first (low 2, K,) the fluid phase 


resistance predominates, and later on (high 2, 


> 


nus 


= 


m 1) the resin phase resistance. 


For extreme values of A, the equilibrium curve 
lies near to the sides of the quadrangle, and can 
be approximated by two right lines ; one from the 
Oy = 1, from there to 


origin to a and one 


‘ From Fig. (14) it is clear that (47) 
and (48) can then be written: 


a ly 1. 


w 


Cy ' Mw , 
kp Sa (50) and . kyS (1 — ax), (51) 


Mt 


since at first 2* — Oand later f(7) — 1 (and y 


vr). 


3. 14, Extremely high equilibrium constant. 


Equations (50) and (51) will give linear relations 


between ¢ and In x or In(1 r). Then. from a 


graph of In a and In(l r) against ¢, one can 
caleulate ky and ky, from the slopes of the two 
lines. 

For normal equilibrium curves, where (50) and 
(51) can not be applied, the main difficulty in the 
interpretation is the fact, that the point («,, y,) 
(or the angle 6) is unknown. In principle the 
inflexion point concentration 1, together with the 


Wevers 


shape of the equilibrium curve could give the 
magnitude of 6, which is suflicient to solve the 
problem. However, it is diflicult to estimate the 
still, 
combining the information provided by the estim- 


place of an inflexion point accurately ; 


ated value of »,, with data from analogous experi- 
ments [e.g. interpreted with (50) and (51)], it is 
possible to make a reasonable estimate of 4. 

Then (4 — 2,) and (y, 
and ky, and k, can be calculated from the slope of 
the break-through curve. 


y) from (45) are known 


6 Kinetic Data From EXxrreRIMENTS 
In Table 1 a summary is given of the systems 


used for kinetic experiments, 


Table 1. 





Ton cachanger Ton pair Anion 
Na’ 
Na 
He 


H 
iH 
Na* 


Amberlite IRC 50 OH 


Dowex 50 » 





Nearly all the experiments were performed in 
“ plastic’ columns; only those with systems 
Na‘ H+ (so,?) Dowex 50 & 8 and 
Ag* — H* (NO,~) —- Dowex 50 « 8 were done in 
glass columns. 

The equilibrium curves of the systems used with 
Dowex 50 « 8 were taken from the first part of 
this paper. The dissociation constants of HOH, 
H,BO, and H,SO, were so far different from that 
of the — COOH groups in the resin that the three 
Amberlite IRC 50 systems could be interpreted 
with relations (50) and (51). 

From the values found for ky a diffusion co- 
eflicient in the resin phase By was calculated with 


kyS d? - ky (l — e)d 


D 
ad 60 10 


(52) 


It has been argued already that the rate relation 
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(38) cannot describe the real 
accurately, That is the reason why a number of 
authors has used other constants in (52); the 
following numerical values are used: 60 or 47? 
[7, 10], 4x? or 48 [6] and 30 [13]. The choice is 
rather arbitrary and not very important, since 
from the simplified description [by (38)] one should 
not expect more than an order of magnitude of B »,. 


Criterion (39) was sufliciently fulfilled for all 


phenomena 


kinetic experiments to be described. 


Kinetic experiments with Amberlite IRC 50 
Table 2 gives a summary of the particle diameters 
and specific surfaces used. 


Table 2. 





Mesh no. 
(ASTM) 


O42 


Effective diameter S.10-2 (n-) 
(mim) 


Unscreened | 89 (indicatedwith * in 
following tables) 
30 0-64 38 
0-42 Bt) 





The results of the experiments are given in 
Tables 3, 4 and 5 for the resin phase resistance for 
H,BO,-, OH 
6 and 7 for the liquid phase resistance for SO, 
and H,BO,~ anions. 


and SO, anions and in Tables 


Table 3. 
(H,BO,~) 


Results of experiments with Na*—H* 
Amberlite IRC 50 (resin phase); Co 
0-134 — 0-137 kg eq./m*. 





S.10>2 
(mot) 


Experiment — Column diameter DBD, lo” 


no. (mm) (m2? /s 
sv 
sv 
sy 
sv 
su 
su 
sv 
sv 
58 
58 


Mean 





Table 4. 
Na* 


Results of kinetic experiments with 
H*(OH~-) — Amberlite IRC 50 (resin phase) ; 
Co 0-31 


0-51 kg eq./m?* 





Experiment Column diameter S.10>- 2 DB, 10 


no. (mm) (m-) (m2? /sec) 


5-03 
5-61 
4-38 
3-81 
471 
3-44 


29° 10 SY 
33 10 xo 
SS 10 SY 
OS sv 
7 sv 
101 8Y 


Mean 





Table 5. Results 
H* —~Na* (SO, ) 


of kinetic experiments with 
Amberlite IRC 50 (resin phase). 





Column 
S.10-2 


no. (mm) (m4) 


Experiment diameter D, 10” Co 


(m* sec) | (kg eq. /m*) 
20-1 
38-5 
37-9 
42°55 
20-8 
38-7 
42-1 





Table 6. 
H+ — Na* 
phase ; Co 


Results of kinetic experiments with 
(SO, ) IRC 50 (liquid 
ca. 0-4 kg eq. m®, S — 89. 10®m~!), 


Amberlite 





ky.10® 


m sec 


. 2 
Ug. 10 


mi sec 


Sh Re 


773 
6-64 
13°31 
5-95 
7-31 
551 


3-45 


0-087 
0-075 
0-150 
0-067 
0-082 
0-2 
0-039 





Not all experimental break-through curves could 
be used for the determination of both k, and kx. 
That depended on the number of observation 
points either in the beginning (for k,) or the end 
(for ky) of the break-through curve. 
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Table 7. Results 
Nat —H* (H,BO, ) 
phase ; Co = 0-134 


of kinetic experiments with 
Amberlite IRC 50 (liquid 
0-137 kg eq. m', S = 89. 108 
m-~') 





ky,.108 U,-104 Re 

176 
135 
OoT2 
0-082 


15-70 
11-08 
6-36 


725 





the numbers of 


defined accord- 


In Tables 6 and 7 values for 
Reynolds and Sherwood are given, 
ing to: 

hy 


Re - 0 
SB, 


S — 
A Sv 


Here D, diffusion coetlicient in the liquid phase 
(m?* sec), 
kimematic the 


viscosity ot licpuacd 


(m* sec), 


For v the value of 10-® m? sec was always taken 
(= viscosity of water at 20 C), because all experi- 
ments were performed at around 20 C. Less simple 
is the choice of a value for B, since different ions 
are diffusing in opposite directions. In general the 
smallest diffusion coeflicient of the two substances 
present was chosen. 

For the Na*—H*(H,BO,~) system the value 
for boric acid was taken, that for borax being 
unknown. Table 8 gives the values used (L.C.T.) 


Table 8. 





System DB, 10° (m? sec) 


thas 


4 
H* (H,BO,~) 
Hit (CI-) 

H* (NO,~) 


at —~ H* (SO2-) | 





In Fig. 15 a double logarithmic graph ts given, 
representing the experimental results of Tables 
6 and 7. The drawn lines in this Fig. can be 
represented by : 


Sh = C Re" 
with n «= 0-8. 


] 
E 
} 


Vin 


“e 
Fic, 15. Sh as a function of Re for the system Na*— H’ 
Amberlite IRC 50. 


Kinetic experiments with Dower 50 « 8 

In view of these results it is assumed for the 
Na*—H*(Cl-) experiments with Dowex 50 « 8, 
that could not be interpreted with (50) and (51), 
that ky and tgd — Cy ky Gok, varied proportional 
tovg’*. Further it was apparent from experiments 
Table 9, that the break-through 
curve at a medium flow rate was practically 


summarized in 


symmetrical. Therefore 6 was taken to be 45° for 
experiment number 31. 

Then it is possible to determine k, and ky from 
the break-through curves and the equilibrium 
curve given in Fig, 3. 

The results are given in Table 9. 

From the mean value of k, follows with (52) 
BD, 
with the value 1-6.10°% m? see of GILLILAND and 
Bappour [6] 8:10-"! — 1-10° m?*/see of 
Boyp and So.pano [2]. 

From two other experimental break-through 


1-2.10-2° m? ‘sec, comparing reasonably well 


and 


curves mass transfer coctlicients were determined ; 
the systems involved were Nat—H* (SO,*>) 

Dowex 50 * 8 Ag* — H* (NO,~) —- Dowex 
50 x 8. The inflexion point concentrations «, 
combined with the equilibrium curves (Figs. 4 and 


and 


5) show that in the first mentioned case the liquid 
side resistance predominates (6 ~ 0°) while in the 
second case 6 = 45°. The slope of the break- 
through curves was determined graphically for a 
few values of x. The results are given in Tables 


10 and 11. 
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Table 9. 


Results of experiments with the system Nat —H*(Cl-) 


Dowex 50 « 8 





Experiment 
no, (m /sec) (see~1) 


5 


21-1 8-06 
12-8 
58 300 
22-6 6-80 
28-1 


oO 
-_- 


625 


—_— = 


26-8 
34-3 


oo Se 


“os 


(m/sec) 


| 
Up 104 | (~) 108 kp. ky, .10° Re 
Mt) 5 


(m/sec) 


O-155 
0-094 
0-042 
0-166 
0-206 
O-O16 
0-197 


324 
2-05 
0-99 
2-88 
15 
0-56 
4-08 
4-97 


5-28 
334 
1-62 
1-69 
“77 
O-o2 
O57 


S-1l 


252 





Table 10. Results of kinetic experiment with the 
Nat —H* (SO?) — Dower 50 & 8. S 
99.10? m-!; glass column, diameter 
U9 1:3.10°3m sec; Re O-13L; Cy 


. . 3 
kg eq. m'*, 


Gj) 


(sec *) 


system 
14mm; 
0-302 





(ky .10® 


(mi sec) 





Table 11. Results of kinetic experiment with the 
Ag’ H (NO, ) Dower 50 & 8 OS 
192.102 m~!; glass column, cross-section 


Vg = #78.10° 4m sec; Re — 00-0249; Cy 


system 
9 

lem-*; 

O-202 


kg eq. m?*. 





| (dw /dt) 
r 13 


| (see*) 


| kp. 108 
(m sec) 


kp. Ww® 
(msec) 


O40) 3-91 
0-50 4-85 
0-75 | 3-64 


9 
1-7.107 Um? see 


Dr 


| 
| 
| 





7. Discussion or tur ResuLtTS 
Resin phase resistance 


From Tables 3, 4 and 5 it is seen that the resin 
phase diffusion coetlicients By, differ rather much 


for the three systems studied, that for borax 
being the lowest, and sulphuric acid the highest. 
A qualitative explanation of this phenomenon can 
be given as follows. 

Upon exchange of Na* and H* ions with SO, 
as the anion, only the small mobile cations are 
The 


taneously with the cations, the salt Na,SO, being 


diffusing. anions need not diffuse simul- 
totally and the acid H,SO, practically dissociated. 
If, however, OH 
HOH, is very weakly dissociated. Consequently, 
the OH 
the Na 

The same applies to the borax experiments 


Na’ 


currently with practically undissociated boric acid ; 


is the anion, the acid formed, 


ions must enter the resin together with 
ions against a stream of undissociated H,O. 
where and H,BO,~ ions diffuse counter- 
only in this case the ions involved are bigger and 
less mobile than with OH>~ anions, resulting in a 
lower value for Dp. 

Conway et al. [3] found a value of 0-4.10-” 
m? see for the Nat —- H* (OH-) — Amberlite IRC 
50 system ; that is an order of magnitude smaller 
They calculated B, 
from batch experiments in which they measured 


than the value from Table 4. 


the progress of this exchange in a stirred vessel. 
In their interpretation they supposed that no 
liquid side resistance was present. However, the 
value of m for this exchange will be extremely high 
in the beginning of the experiment ; therefore the 
liquid phase resistance is bound to predominate 
in the first stages of the experiment. That means 
that the results of Conway et al. are not reliable 
since the specific properties of the resin make it 
unsuited for their method of measuring. 
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Liquid phase resistance 

In Fig. 16 all liquid side kinetic measurements 
are summarized in one graph, giving Sh as a 
function of Pé (Pé — vg SB,). 

The following points 
Fig. 
(a) The points do not fall on one line as was to be 
An for this 
spread is probably the fact that the values 


are evident from this 


expected. important reason 
taken for B, were wrong. The same reasoning 
explaining the difference in Dy applies here 
too. As a consequence the diffusion coeflicient 
for the borax system—and to a smaller extent 
for the sulphuric acid system too— should be 
much smaller relative to that for the chloride 
The effect of this is that the pomts 


in Fig. 16 are drawn together considerably. 


system. 


(b) Notwithstanding the spread in the results, it 
is clear from the 
of the Reynolds, 
cidedly higher than 0-5 [14, 19] and more 


graph that the exponent » 
or Péclet, numbers, is de- 


nearly 0-8. The latter value agrees with that 


found by other authors, who performed 


kinetic experiments with ion exchangers [6, 
18]. 


== 








—_ 


Pé 


Fig; 16, Sh as a function of Pé, 


OF THE EXPERIMENTAL 
TECHNIQUE 


8 DrTAILS 


Theexperimental procedure for measuring break- 
through curves is simple and will not be given here 


in much detail. Only the most important points 
will be mentioned in the following section. 

The equilibrium experiments described were 
performed with the cation exchanger Dowex 
50 « 8. The kinetic experiments were performed 
with the same resin and also with Amberlite IRC 
50. 
sieved and a certain sieve-fraction was used, 


In most cases the exchanger particles were 
The 
particle diameter is not very important for the 
equilibrium experiments but it is important for 
the kinetic experiments. It is necessary to discuss 
the consequences of the swelling and shrinkage 
occurring in a bed of ion exchanger particles, 
These volume changes are much more pronounced 
in the case of Amberlite IRC 50 than for Dowex 
HO) & 8. 

As a rule the exchange resin shows a shrinkage 
when there occurs an exchange, in which the 
displacing ions have a greater affinity to the resin 
than the displaced ions. Consequently a bed of ion 
exchange particles will normally shrink if K, > 1 
(sharp front) and swell if K,- 
break-through). 
important consequences for the measurement of 


1 (“ equilibrium ” 
This behaviour has some very 


break-through curves. 

If in a percolation experiment A, > 1, a sharp 
break-through curve will result and shrinkage will 
occur. If the packing of the bed is not absolutely 
regular, the flow rate will be somewhat greater in 
the places with high porosity and smaller where the 
porosity is less. So, at the places with high 
porosity, the front will run ahead, causing more 
shrinkage. 

In this way small and often inevitable differ- 
ences in the bed porosity will grow more pro- 
nounced during the experiment, which means that 
a flow as deseribed here is in principle unstable. 
Therefore experiments with sharp break-through 
curves should be mistrusted principally if no 
attention has been paid to this aspect. 

If, on the other hand, K, <1, the same reason- 
ing leads to the conclusion that the results of a 
reverse or “ equilibrium ” break-through curve are 
much more to be trusted. Further the longitud- 
inal mixing, resulting from an unequal flow rate, 
is much more serious in the case of the sharp 
break-through curve with great concentration 
gradients in the direction of flow than for the 
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equilibrium break-through where the concentra- 
tion changes only slowly. 

Therefore the experiments described in the first 
part of the article were mostly performed with 
normal glass columns, but the sharp break-through 
curves used in the second part, describing the 
kinetic results, were measured using an other type 
of column. 

The wall A of this column (Fig. 17) is made of 
slightly extensible * plastic ” material (e.g. silicone 
rubber). The lower part of this column is clamped 
(clamp D) on a short piece of glass tubing C, 
bearing a fritted 
ending ina capillary. The construction is apparent 
The bed of exchanger particles in 


glass disk B and a glass cock 


from Fig. 17. 
A can expand laterally producing a certain stress 
in the column wall. This stress will, on subsequent 
shrinkage, prevent any uneven packing in the bed 
practically completely. The break-through curves 
measured with this type of column were far better 


reproducible than those taken with glass columns. 


For a more complete description of these phenom- 
ena reference is made to [20]. 


Fic, 17, “ Plastic ” column. 


It is 
through curves with a swelling bed, thus prevent- 
ing the deleterious effects of shrinkage on a sharp 
break-through curve. case it is 


also possible to measure sharp break- 


For such a 


necessary to make use of the weakly acidic 
character of the exchange groups of e.g. Amberlite 
IRC 50, 


dissociated if they are in the H* form. 


These carboxylic groups are weakly 
Therefore 
during the exchange of e.g. Na*, against H*, the 
following reaction occurs. 


R COOH = R COO 
H+, + Nat, —— Nat, 


H*» 
H*, 


(34) 
(35) 


The dissociation reaction (34) has such a small 
equilibrium constant A, that the H> ions are far 
more strongly bound to the resin than Na* ions. 
Therefore it is impossible to displace the H*, ions 
from the resin by a solution of NaCl or Na,SO,. 
It is possible, however, to draw reaction (35) to 
the right hand side by using an anion of a very 
weak acid. Then the following reaction will occur 
in the solution phase : 


H+, - 


If the 
smaller than A 


dissociation constant A, of (36) is much 

; of (34) the equilibrium of (35) will 
lie at the right hand side (or A; > 1). This situ- 
ation will occur if for Z the anion H,BO, or 


OH 


But since the hydrogen in the resin phase occurs 


is chosen. 


almost completely in the undissociated form, so 
without water of hydration the resin will swell 
considerably when the hydrogen is displaced by 
the hydrated Na* ions. So in this case a sharp 
break-through curve (A; 
a swelling bed, 


- 1) will be attended by 
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NOTATION 


absolute concentration of one ion species in 
the fluid (kg eq. m*) 
value of C at the interface between two 
phases (kg eq. /m*) 

total concentration in the fluid (kg aq./m?) 
constants in (41) 

constant 

particle diameter (m) 
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= diffusion coefficient in the liquid phase 
(m?®, sec) 
= diffusion 
(m?® sec) 


coefficient in’ the resin phase 

= function of a in (5) 

- area of the cross-section of the column (m*) 
reaction rate constant in (37) (m* kg eq. 

sec) 

partial mass transfer coefficient in the liquid 

phase (m_ sec) 

partial mass transfer coeflicient in the resin 
phase (msec) 

dissociation constant in (20) (kg mol, m®*) 

dissociation constant of the active groups 
in the resin (kg mol m*) 

total mass transter coefficient for the liquid 
phase (msec) 

total mass transfer coeflicient for the resin 
phase (m/sec) 


equilibrium constant for uni-univalent ex- 


Rp 


change in (13) 


equilibrium constant for uni-divalent ex- 


change in (33) 


= dimensionless measure of the volume of 
fluid that has flowed through the column 


special value of 1, where / 1/|4e (1 — «)) 


} 


slope of the line connecting the points 
(0, O) and (a2, 4) 


= exponent of the Reynolds number 


absolute concentration of one ion species in 


the resin (kg eq. m® bed volume) 


value of q at the interface of two phases 


(kg eq./m® bed volume) 


total capacity of the exchanger (kg eq. /m* 


bed volume) 


volumetric flow rate of the fluid (m®* sec) 


partial resistance against mass transfer in 


the liquid phase (m*sec /kg eq.) 


(1) 
2) 
[3] 
[4] 
(3) 
(6) 
(7) 
[8] 
[9] 
{10} 
(11) 
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Shorter Communications 


The diffusivity of chlorine in water 


( Received 17 November 1958) 


We hive recently determined the diffusion coefticient for 
chlorine in water by measuring the rate of absorption of 
pure chlorine gas into water in a wetted wall column, The 
measurements were performed at atmospheric pressure and 
at temperatures ranging from 10 to 35 ©. 

The experimental technique was exactly the same as 
employed by Nusstne ef al. (1) and the hydrodynamic con- 
ditions of the laminar film were accurately known, As has 
been explained in [1] and in many other recent papers on 
the same subject, the rate of absorption can be calculated on 
the basis of the * penetration ~ theory, From this it follows 
that for physical absorption the absorption rate is propor- 
tional to c* , D, the proportionality factor being fixed by 
the width and the height of the falling film and its surface 
velocity. 

In this investigation the theory of physical absorption 
has been applied to the absorption of CL, into water, 
although it is known that under equilibrium conditions the 
chlorine in water is partly hydrolysed. Thus we find from 
the data collected by Apams and EnoMonps [2] that of the 
total amount of chlorine dissolved about 78 per cent is 
unhydrolysed at 10 C, 71 per cent at 20 C and 63 per cent 
at 30°C. However, from «a paper by Sumov and 
SOLODUSHENKOV [3] it appears that at room temperature 
the hydrolysis reaction is completed well within 0-1 see. 
Since the contact times in our film absorber were between 
0-3 and 1 sec, we have assumed chemical reaction equil- 
ibrium throughout the liquid phase. In that case, the 
theory of physical absorption may be applied with c* = the 
equilibrium solubility of chlorine in water (both hydrolysed 
and unhydrolysed) and D = the diffusivity of the equil- 
ibrium mixture of Cl,, HOC] and HCI in water. 

In Table 1 the values obtained for c*, D are given 
which have been corrected to a chlorine pressure of 760 mm 
5 experimental runs 


Hg. Kach value is the average of 3 


with different liquid velocities and in each run the film . 


height was varied by a factor of 6. The maximum deviation 
of the individual c* , D values for one temperature from 
the average values indicated in Table 1 did not exceed 
% percent. The mean deviation is about + 1-5 per cent. 
For calculating the values of D at the different tem- 
peratures, the solubility data of Apams and EpMoOnps [2] 
have been used. These are given in Table 1 together with 
the calculated values of D. If the solubilities between 
10 and 25°C from the paper by Warrney and Vivian [4] 
are employed, the value of D becomes 5 per cent smaller 
at 10°C and 5 per cent larger at 25°C. It appears that the 
reliability of the diffusivities obtained in this investigation 
depends mainly on the solubility data adhered to. 


Table 1. 


Experimental values of c* yD and derived values 


of D 





c* (2) + D 
(kg m*®) | (10-9 m? sec) 


ety D t 
(10-4 kg m?® see?) | 


Temp. 
() 


ool 
O08 
20 


OS 9-76 
S4 9-05 
70 7-76 
is 7-37 


1: 
1- 
“0 1-32 
1- 
1- 
1- 


1-0 : 
1s 
18-3 
20-0 
224 
25 

30 


22 
a) 
“i 654 
5-86 
5435 


$2 
2 


SI 


= 


* * * 
meee nw wwe a 


- 
— 





+ At 760 mm Hy pressure. 


The diffusivity values of Table 1 have been plotted 
against temperature in Fig. 1, together with the data 
obtained by Davipson and CuLLeN [5] from absorption 
experiments on wetted spheres, As to the dependence on 
temperature, their diffusivities for chlorine were rather 
discordant with those obtained for other sparingly soluble 
non-reacting gases with the same method. This is not the 
case with the chlorine diffusivities presented in this com- 


2.0 , , — a 


"a 





——1= 0(107*m? sec-') 
- 
wr 





a 
o 





poper 











i 
10 20 30 40 


———-= temperature (°C) 











Diffusion coeflicient of chlorine in water; values 


from Table 1 vs. temperature. 
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munication, It appears that the Stokes — Einstein relation- — Laboratorium voor Fysische Technologie H. Kramers 
ship holds well for these data ; the average value of Du T = Technische Hogeschool R. A. Dovueias 
is (423 4 0-07) lO N°K-!, Delft, Netherlands R. M. ULMANN 
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Book Reviews 


The Kinetics of Vinyl Polymerization by Radical 
Mechanism. C. H. Bamrorp, W. G. Bars, A. D. Jenkins 
and P. F. Crvax. Butterworths, London, 1958. 318 pp. 


Ss.  SS.80. 


Free radical polymerization has been the subject of 


extensive scientific research for a period of somewhat 


more than 20 years. 
phenomena involved has been established, together with a 


During this time the nature of the 


detailed mathematical treatment of the several processes 
which are believed to take place. Unlike the other pro- 
cesses for making high polymers the free radical process 
is now capable of quantitative description through the 
absolute velocity constants for the different steps involved. 

The time is therefore ripe for the production of a single 
matter of 
satisfaction for all that 
Dr. Bamrorp and some of his colleagues have found the 
The work of the Maidenhead laboratory 


work on this subject, and it can only be a 


workers in the field to see 


time to write it. 
has made a number of important contributions to the 
understanding of free radical processes, especially in 
connection with the problems of trapped radicals and 
absolute velocity constants, and these are naturally well 
represented, (though not over represented) in the text. 

The book covers all the important fields of radical 
polymerization, It begins with techniques of investiga- 


tion and then deals with the measurement of absolute 


velocity constants for different monomers. Subsequent 
chapters show how these quantities work out in terms of 
and heterogeneous 


homogeneous polymerizations, 


Evidence is presented showing how the reactions are 
affected by the of free 


particles or by their fixation in gels or precipitated polymer 


isolation radicals in emulsion 


particles. Other chapters deal with copolymerization, 
of which a useful condensed presentation ts given, and 
the kinetics of inhibition and gelation. 

\ book with this title naturally contains a good deal 
of mathematics, but it may also be read with advantage 
by those who do not wish to go through all the mathe- 
matical arguments in detail. The inclusion of an author 
index would have been an advantage, as the book should 
become an important reference source for future workers 
in the field, 


Kh. N. Hawarp 


Kurzes Handbuch der Brennstoff- und Feuerungs- 
technik. W. Guz: Springer Verlag, Berlin Gottingen 
Heidelberg, 1953. (2. Aufl.) 583 8. DM 45 


Das 1942 in erster Auflage erschienene und seit 1952 in 
2. Auftlage vorliegende Buch gibt auf knappen Raum 


(575 5.) in 4 Hauptabschnitten 
gesetze, die Brennstoffe, Verbrennung, Vergasung, Ver- 


Physikalische Grund- 
brennungs- und Vergasungsvorgiinge eine umfassende 
Einfiihrung in alle Fragen und Probleme, die bei der 
praktischen Verbrennung der herkOmmilichen gasfOrmigen, 
fliissiven und festen Brennstoffe auftreten. Die moderne 
Feuerungstechnik wird von dem Verfasser als angewandte 
Physikalische Chemie aufgefasst, die neben der Thermo- 
dynamik in immer steigendem Masse genaue Kenntnis 
der StrOmungslehre verlangt. 

Finftel 
seines Buches den diesbeziiglichen physikalischen und 


In diesem Sinne widmet der Autor etwa ein 


physikalisch-chemischen Grunderscheinungen. In einer 


konzentrierten Aufzihlung aller wichtigen Gesetzmiissig- 
keiten wird dem Feuerungstechniker gezeigt, was es evtl. 
alles zu beriicksichtigen gilt. In cinem Kapitel iber die 
Brennstoffe tindet sich cine Beschreibung der physikalischen 
und chemischen Stoffeigenschaften sowie der Bewertungs- 
moOglichkeit durch Verbrennungswiarme, Reaktionsfihigkeit 
und Ziindeivgenschaften. Auf 60 Seiten sind ausserdem 
im einzelnen die Aufarbeitungs (und Veredelungs) methoden 
dargestellt. Im Mittelpunkt des Buches steht als besonders 
wert voller Beit rag cine cingehende und kritische Einfiihrung 
Berechnungsmoglichkeiten des 


in die verschiedenen 


Verbrennungsvorganges unter besonderer Berticksichti- 
gung des I-t Diagramms, zu dessen Entwicklung der 
Verfasser selbst entscheidend beigetragen hat. 
Bemerkenswert ist die Fille der zusammengetragenen 
Zitate, die diirften, da der 
Verfasser durch seinen langjahrigen Aufenthalt in 
USA guten Einblick in die dortige auch weniger leicht 
zuginyliche | Literatur hatte Etwa 200 treffend aus- 


besonders vollzihlig sein 


gewihlte Abbildungen und rund 100 Zahlentafeln unter- 
stiitzen die konzentrierte Darstellung. 

Man darf den Versuch als gelungen betrachten, durch den 
hier dem Praktiker ein Buch an die Hand gegeben wird, 
das bis zu den Grundlagen fiihrt und nicht in einer Ausseren 
steckenbleibt. Die vertiefte Darstellung 
bietet sicherlich auch dem Forscher manche Anreguny. 


Beschreibuny 


Es wiire zu wiinschen dass es dem Verfasser in Fort- 
fihrung seiner verdienstvollen Arbeit gelingen moéyge, das 
Buch stindig auf neuzeitlichem Stand zu halten. Bei 
einer Neuauflage wiire dann evtl. die chemische Kinetik 
mehr als bisher zu beriicksichtigen. Das immer stirker 
werdende Eindringen kinetischer Betrachtungen in die 
praktischen Anwendungen lisst dieses Gebiet zu einer 
eigenstandigen 
mungslehre ete. 
tiberlegt werden, ob ein Abschnitt fiber die wichtigsten 


Disziplin neben Thermodynamik, Stré- 


werden. Weiterhin koOnnte vielleicht 


Grundergebnissne der Flammenforschung mit in das 
Kapitel iiber die Grundgesetze aufzunehmen wiire. 
F. Frerrine. 
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In Fig. 1 on p, 82 x = O should read vy = 0 





